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Millimeter  Wave  Scatter  and  Attenuation 
Measurements  on  Snow  Slabs 


1.  INTRODUCTION 

Millimeter  wave  component  technology  has  advanced  to  a  state  where  a  large 
number  of  military  system  applications  are  coming  under  consideration.  The  ad¬ 
vantages  of  small  size,  high  antenna -beam  resolution,  and  large  bandwidth  are 
widely  recognized,  particularly  onboard  space -limited  aircraft  and  missiles.  Propa¬ 
gation  phenomena  such  as  high  attenuation  in  the  molecular  resonance  bands  and  in 
heavv  precipitation  are  not  a  problem  in  many  tactical  applications.  In  fact,  covert¬ 
ness  beyond  the  intended  operational  range  is  considered  an  advantage  in  many 
cases.  Imaging,  ranging,  and  homing  millimeter  wave  sensors  interact  with  natural 
backgrounds.  It  has  been  found  that  old  granular  snow  covering  the  ground  can  have 
unusually  strong  reflective  properties  at  millimeter  wavelengths.  The  implications 
of  a  high  reflectivity  of  metamorphic  snow  must  be  of  concern  to  military -system 
designers. 

A  recent  survey  by  Suits  and  Guenther^  on  the  scattering  of  natural  backgrounds 
in  the  millimeter  wave  region  showed  only  limited  data  in  the  present  literature. 
Measurements  at  shorter  wavelengths,  on  old  granular  snow  in  particular,  arc 
scarce  or  non-existent.  A  further  drawback  on  snow  measurements  quite  often  is 


(Received  for  publication  17  September  1981) 

I.  Suits,  G.  H.  ,  and  Guenther,  B.  D.  ( 1978)  Targets  and  Backgrounds,  U.S.  Army 
MIRADCOM.  Tech.  Rep.  T -78-7  1. 


the  lack  of  an  adequate  description  of  the  meteorological  conditions  and  snow  com¬ 
position  under  which  the  scatter  data  were  obtained. 

•> 

l'he  University  of  Kansas-  lias  reported  on  snow  backscatter  measurements, 
conducted  in  the  winter  of  1975  over  the  frequency  range  1  to  18  GHz.  During 
February  and  March  of  1977  further  measurements'*'^  of  snow  backscatter  in 
Colorado  extended  the  frequency  range  to  include  1  to  18  GHz  and  35  GHz.  Data 
were  acquired  at  numerous  frequencies,  polarizations,  and  grazing  angles  for  a 
variety  of  carefully  observed  snow  conditions.  F lie rbruch  et  al3  at  the  National 
Bureau  of  Standards  have  used  frequencies  between  0.  5  and  18  GHz  to  probe  the 
stratigraphy  of  snow  based  on  backscatter  techniques.  Antennas  mounted  directly 
above  the  snow  detected  reflections  from  depths  up  to  4  m  depending  on  frequency. 
Hoekstra  and  Spanogle  determined  the  scattering  from  snow  in  New  Hampshire 
at  low  grazing  angles  for  10  and  35  Cl  Hz.  Prompted  bv  earlier  Honeywell  findings 

of  snow  effects  on  experimental  millimeter  wave  sensors,  the  Georgia  Institute  of 

7  8 

Technology  conducted  measurements  '  at  35  and  95  GHz,  again  in  Colorado  during 
the  later  winter  of  197  5/7  6.  Despite  limitations  in  the  experiment,  evidence  was 
found  of  some  of  the  peculiarities  of  metamorphie  snow  at  millimeter  wavelengths: 
its  ability  to  produce  strong  scattering  in  the  frozen  state,  the  existence  of  changes 
in  reflectivity  associated  with  the  diurnal  melt-freeze  cycle,  and  the  high  attenua¬ 
tion  of  signals  passing  through  snow,  suggesting  that  scattering  is  a  near -surface 

9 

effect.  Ulaby  et  al  have  published  a  state-of-the-art  review  of  the  microwave 
backscatter  and  emission  behavior  of  snow.  Trebits  et  al10  in  another  review  of 
millimeter  wave  background  scatter  compare  snow  data  obtained  by  the  University 
of  Kansas  and  the  Georgia  Institute  of  Technology.  They  point  out  that  a  substantial 
increase  in  snow  reflectivity  has  been  observed  when  going  from  10  to  35  GHz. 

In  order  to  better  understand  these  phenomena,  and  to  obtain  quantitative  data 
at  35,  98,  and  140  GHz,  the  Propagation  Branch  at  RADC  has  conducted  a  measure¬ 
ment  program  during  the  winter  of  1977/7  8  and  investigated  backscatter,  bistatic 
scatter,  and  attenuation  characteristics  of  metamorphie  snow  at  these  frequencies. 


2.  SCOPE  OF  MEASUREMENT  PROGRAM 

A  choice  had  to  be  made  with  regard  to  the  frequencies  at  which  to  obtain  data. 
Component  technology  and  hence  exploratory  system  development  is  most  advanced 
in  the  35-GHz  region  and  to  a  slightly  lesser  extent  in  the  94-GHz  region.  The 
status  of  140  GHz  and  higher -frequency  technology  must  still  be  considered 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  73.  ) 
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experimental.  Based  on  the  current  need  to  establish  the  tradeoffs  between  35-GHz 
and  94-GHz  systems,  and  a  possible  requirement  for  140-GHz  data  in  the  near 
future,  the  determination  was  made  to  derive  comparative  data  at  these  three  fre¬ 
quencies. 

In  many  of  the  prospective  applications  of  millimeter  waves,  imaging  of  natural 
backgrounds  or  of  man-made  targets  against  natural  backgrounds  plays  an  important 
role.  Imaging  or  target-to-background  studies  involve  complexities  that  would 
distract  from  a  study  of  basic  millimeter-wave/snow  interaction.  No  effort  was 
undertaken  during  the  present  phase  of  the  program  other  than  the  direct  measure¬ 
ment  of  snow  scatter  and  attenuation  parameters. 

2.1  Hie  Sraltrronn'ter-Traiisinissoinetpr 

During  the  planning  stages  of  the  experiment  it  became  obvious  that  suitable 
commercially  developed  equipment  was  unavailable  to  ItALK’.  The  decision  was 
made  to  design  and  construct  short-range  scatterometers  in-house  which  would 
permit  the  measurement  of  baekscatter  coefficients  and  bistatic  scatter  coeffi¬ 
cients  as  well  as  attenuation  coefficients  on  samples  of  metamorphic  snow, 
removed  from  natural  snow  surfaces.  In  later  stages  of  metamorphism,  this 
can  be  accomplished  without  disturbing  the  snow  in  a  noticeable  way.  The 
main  objectives  in  the  equipment  design  were:  (a)  flexibility  in  measuring 
a  wide  range  of  propagation  parameters  on  snow,  (t>>  ruggedness  for  operation 
in  the  field,  and  (c)  the  use  of  components  available  in-house  to  accelerate  equip¬ 
ment  construction  in  order  to  perform  measurements  during  the  winter  1977  / VH. 

Standard  gain  horn  antennas  for  transmission  and  reception  at  the  three  fre¬ 
quencies  permit  fat— field  measurements  at  ranges  of  1  m  with  3-dB  beam  diameters 
of  approximately  0.  24  m.  This,  combined  with  the  fact  that  old  granular  snow  has 
sufficient  rigidity  to  permit  the  removal  of  large  frozen  slabs  from  the  ground,  led 
to  the  development  of  a  sample  measurement  technique  incorporated  in  the  equip¬ 
ment  shown  in  Figure  1.  A  sled-like  platform  was  constructed  from  plywood  with 
metal-underlaid  runners  to  permit  easy  movement  on  snow.  A  double  wooden  boom 
is  rigidly  attached  to  the  platform,  extending  upward  at  45  deg.  At  its  upper  end, 
this  boom  carries  a  wooden  basket  to  hold  a  snow  sample  of  0.  76  m  length  and 
0.56  in  width.  The  sample  holder  rotates  around  a  horizontal  axis.  A  double  wooden 
boom  extension  pivots  around  the  same  axis.  On  the  end  erossmember  of  the  boom 
extension  the  receivers  are  mounted.  The  transmitters  are  located  on  a  eross¬ 
member  at  the  lower  end  of  the  rigid  boom  near  the  sled  platform.  Both  the  trans¬ 
mitting  and  receiving  antennas  are  pointing  at  the  common  axis  of  rotation  and  hence 
at  the  center  of  the  snow  sample.  This  is  independent  of  the  angular  position  of  the 
snow  surface  relative  to  the  transmitter  and  receiver  antenna  directions.  In 


the  configuration  shown,  the  transmitter-target-receiver  geometry  approaches  that 
for  backscatter.  The  maximum  angle  between  transmitted  and  scattered  signal 
directions  is  WO  deg  with  transmitting  and  receiving  antennas  almost  touching  each 
other.  The  snow  sample  in  the  basket  is  inverted  from  its  natural  position  and  is 
prevented  from  falling  out  by  a  woven  net  of  waxed  string  or  nylon  fishing  line.  A 
major  advantage  of  this  design  results  from  the  fact  that  for  the  backscatter  and 
most  bistatic  scatter  geometries  the  antennas  are  directed  toward  space.  The 
necessary  support  structures  are  made  of  wood,  again  minimizing  unwanted  scatter 
signals.  Other  advantages  of  the  design  include  the  ease  and  repeatability  with 
which  specific  scatter  geometries  can  be  configured  and  also  transmission  measure 
ments  through  snow  be  made. 


Figure  1.  Snow  Scattcrometer,  Sample  Backscatter  -Mode 

The  transmitter  and  receiver  antenna  beams,  pointing  toward  the  axis  of  rota¬ 
tion  of  the  sample  holder,  were  fixed  for  the  duration  of  a  scan,  during  which  the 
snow  sample  was  slowly  rotated  by  a  motor  and  gear  drive.  This  resulted  in  a 


range  of  g razing  angles  on  tin-  snow  slab,  limited  m  tin-  backscattcr  mode  on  one 
end  by  the  length  of  the  antenna  "footprint"  exceeding  the  length  of  the  slab.  On 
the  other  end.  tin'  basket  was  stopped  at  a  near- vertical  orientation,  beyond  which 
thi'  snow  slab  is  in  danger  of  falling  out. 

In  the  bistatie  scatter  mode  (figure  2'.  the  scatter  angle  between  the  trans¬ 
mitted  and  received  beams  is  reduced  hv  swinging  the  upper  portion  of  the  boom 
away  from  its  tucked-in  position  in  Figure  1.  1'he  minimum  usable  scatter  angle 
is  reached  w  lien  the  transmitting  antenna  footprint  matches  ’he  snow  si  ib  in  length 
and  the  receiving  antenna  infercep’s  ■  sign.!  re  radiated  .-p.  ‘  angentiallv  from  *h< 
snow  surface.  1  nder  these  conditions,  the  t  r  u.smitti  r  .ad  rc.-.-n or  be im  ingles, 
measured  n-l.iti-.  e  to  the  sur!  >ee  "1  the  snov. .  .re  uriequ  o.  i:i  order  to  obtain  a 
useful  scan  range  in  the  bistatic  mode,  the  scattering  angle  must  be  larger  than  the 
minimum.  In  general,  it  is  chosen  large  enough  so  that  the  case  of  equal  trans¬ 
mitter  and  receiver  beam  angles  is  included  in  the  scan,  i  rum  an  applications  viev 
point,  unequal  angles  are  possible  and  thervlore  of  interest.  During  a  scan,  the 
angle  between  the  boom  sections  is  constant  while  the  basket  rotates.  This  means 
that  transmitter  ,.nd  receiver  beam  angles  change  in  opposite  directions  with  then- 
sum  (the  scatter  angle)  held  constant. 


/ 


Figure  2.  Snow  Scatterometer.  Sample  Itistatic  Scatter  .Mode 


I  1 


establish  a  flow.  Subsequently,  this  would  lead  to  a  concentration  of  water  near 
the  surface.  Measurements  in  the  sample  mode  were  generally  confined  to  frozen 
metamorphic  snow  slabs,  the  type  that  gives  rise  to  the  most  severe  reflections  on 
the  ground. 

in  Figure  4  a  simplified  diagram  presents  the  basic  transmitter  and  receiver 
building  blocks  and  a  schematic  view  of  antenna  and  basket  arrangement.  There 
are  throe  individual  transmitters  mounted  side  by  side  on  an  adjustable  board  be¬ 
tween  the  members  of  the  boom.  The  transmitters,  as  is  seen  in  Figure  5,  con¬ 
sist  of  microwave  sources  and  horn  antennas.  They  can  be  rotated  about  their 
longitudinal  axes,  to  permit  transmission  at  vertical  and  horizontal  polarization. 


Figure  4.  Scatterometer  Block  Diagram 


Klystron  oscillators  (VC  7  IOC  and  VKT  2123A)  are  used  at  98  and  140  GHz.  A  Gunn 
oscillator  (VSA  9010EW)  serves  as  the  35-GHz  source.  Power  output  at  all  fre¬ 
quencies  is  of  the  order  of  1  mW.  The  klystrons  were  operated  at  the  lowest 
beam  voltage  (1500  V)  at  which  oscillations  were  reliably  sustained.  Sufficient 
dynamic  range  for  snow  slab  measurements  was  obtained  with  reduced  power  output. 
Lower  input  power  improved  both  frequency  stability  and  reliability  of  the  mm-wavo 
tubes.  The  transmitters  are  powered  by  separate  power  supplies  (Figure  4)  and 
are  on/off  modulated  by  a  common  squarewave  modulator  at  a  1-kHz  rate.  The 
transmitters  use  standard  gain  horn  antennas  at  35  and  98  GHz  and  a  scalar  feed 
horn  at  140  GHz.  The  footprint  on  the  snow  slab  at  perpendicular  orientation  (approxi¬ 
mately  0.  24  m  in  diameter  between  the  3-dB  points  of  the  individual  beam  pattern) 
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is  well  contained  within  the  slab  area  of  0.  56  X  0.  i 6  m.  Beamwidths  of  trans¬ 
mitting  and  receiving  antennas  are  listed  in  Table  1. 


Figure  5.  Closeup  View  of  Transmitters 


Table  1.  E-  and  H-Plane  Beamwidths  of  Transmitting  and 
Receiving  Antennas 


F  requency 
(GHz) 

Transmitt 
Calc.  3-d 

ing  Antenna 

B  Beamwidth 

Receivin 
Calc.  3-dF 

g  Antenna 
Beamwidth 

0E  (deg) 

(deg) 

HE  (deg) 

dH  (deg) 

35 

9.  1 

6.  8 

13.  9 

11.  8 

98 

9.  3 

7.  0 

12.  1 

9.  9 

140 

7,  2 

7.  2 

8.  5 

7.  0 

Based  on  microwave  components  available  and  the  goal  to  keep  the  equipment 
as  simple  as  possible,  a  single  heterodyning  receiver  was  designed.  The  three 
frequencies  are  received  in  a  modified  Ka-band  waveguide  mixer  (TRG  A965  with 
Schottky  barrier  diode)  and  downconverted  to  IF  frequency  by  either  using  the  local 
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oscillator  (L.  O,  >  fundamental  frequency  or  its  third  or  fourth  harmonic.  The  choice 
of  98  GHz  rather  than  the  center  frequency  of  that  atmospheric  transmission  window 
is  explained  by  the  need  to  tune  a  single  local  oscillator  to  within  IF  range  of  the 
three  transmitted  frequencies.  The  relative  insensitivity  to  frequency  of  the 
scattering  and  attenuation  effects  observed  with  this  equipment  justifies  the  receiver 
design.  The  specific  Gunn  oscillator  used  here  (Varian  VSA  901 OEW)  was  found  to 
tune  satisfactorily  over  a  frequency  range  from  32.6?  to  35.0  GHz,  substantially 
exceeding  its  specified  tuning  range.  With  100  MHz  the  IF  center  frequency,  the 
following  exact  frequency  combinations  were  used:  (a)  transmit  frequency 
140.  1  GHz,  L.  O.  Frequency  140.  0  GHz  or  fourth  harmonic  of  35.  0  GHz,  (b)  trans¬ 
mit  frequency  98.  1  GHz,  L.  O.  frequency  98.0  GHz  or  third  harmonic  of  32.6?  GHz, 
and  (c)  transmit  frequency  of  35.  1  GHz,  L,.  O.  frequency  35.  0  GHz.  An  intermediate 
frequency  bandwidth  of  20  MHz  generally  accommodates  oscillator  drift  associated 
with  temperature  changes  in  the  outdoor  environment.  A  picture  of  the  receiver 
frontend  is  shown  in  Figure  6.  The  unit  is  mounted  on  an  adjustable  board  at  the 
outer  end  of  the  boom  extension.  As  a  unit,  the  receiver  frontend  can  be  rotated 
by  90  deg  around  its  antenna-mixer  axis  to  receive  signals  of  either  vertical  or 
horizontal  polarization.  It  can  also  be  mounted  in  three  laterally  displaced  posi¬ 
tions  for  proper  lineup  with  each  transmitter  antenna.  A  standard  M-band  gain 
horn  with  an  M/K  -band  transition  as  shown  in  Figure  6  is  used  for  reception  at 
both  98  and  140  GHz.  At  the  higher  frequency  the  waveguide  is  oversized.  The  re¬ 
ceiver  beam  pattern  at  140  GHz  was  measured  and  found  to  agree  with  that  calcu¬ 
lated  in  Table  1.  A  standard  K  -band  gain  horn  is  used  to  receive  the  35-Gliz  signal. 

The  receiver  frontend  consists  of  the  antenna,  mixer,  L.  O. ,  bias  supply,  and 
a  35-dH  IF'  preamplifier.  The  remainder  of  the  receiver,  mounted  on  the  back  of 
the  sled,  consists  of  components  shown  in  the  block  diagram  (Figure  4).  Two  cas¬ 
caded  balanced  mixers  (Z-MATCH  DBM-100B)  are  connected  as  electronically 
controlled  variable  attenuators  with  a  range  in  excess  of  7  0  dB  for  automatic  gain 
control  (AGO.  They  are  followed  by  another  35-dB  IF'  amplifier.  After  detection, 
the  1-kHz  squarewave  modulated  signal  is  amplified  in  an  HP  415E  Standing  Wave 
Indicator  of  15  Hz  bandwidth.  The  rectified  i-kiiz  output  signal  drives  the  AG C 
amplifier,  whose  gain  curve  was  adjusted  for  near- logarithmic  receiver  response 
over  a  50-dB  input  signal  range.  The  rectified  l-kllz  signal  also  drives  the  'i  -axis 
of  a  Moseley  XV -recorder.  The  X-axis  signal  is  a  potentiometer-derived  dc  voltage 
proportional  to  the  angular  position  of  the  snow  basket.  The  time  constant  of  the 
receiver  is  essentially  determined  bv  the  mechanical  time  constant  of  the  \1  -  re¬ 
corder.  The  recording  speed,  that  is,  the  speed  of  rotation  of  the  equipment  was 
empirically  set  so  that  the  recording  system  could  trace  the  actual  signal  excursions. 
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Figure  6.  C’loseup  View  of  Single  Receiver  Krontend 


2.2  System  Calibration 

Attempts  were  made  to  refer  the  snow  scatter  signals  to  the  known  return  from 
corner  reflectors  or  small  metal  plates.  The  results  were  generally  not  as  satis¬ 
factory  as  those  obtained  with  a  flat  metal  plate  large  enough  to  exceed  the  size  of 
the  antenna  beam  cross  section.  The  similarity  of  this  plate  to  an  infinite  plane 
and  the  image-source  concept  were  investigated  in  some  detail.  Signal  levels  were 
compared  between  the  equipment  configuration  in  Figure  3  with  the  sample  holder 
empty  and  that  of  Figure  1  with  a  flat  metal  plate  of  0.  56  X  0.  76  m  filling  the  sample 
holder  and  the  boom  extension  swung  into  the  170-deg  scatter  angle  direction.  With 
the  plate  oriented  perpendicular  to  the  line  dividing  the  angle  between  the  trans¬ 
mitter  and  receiver  directions,  signal  levels  during  system  checkout  agreed  to 
within  0.  5  dB  in  both  cases.  This  was  found  at  all  frequencies. 

A  schematic  view  of  the  image-source  calibration  procedure  is  presented  in 
Figure  7.  Two  measurements  are  performed  in  sequence.  First,  the  return  from 
the  "infinite"  metal  plate  at  range  r^  from  the  scatterometer  is  determined.  For 
the  purpose  of  illustration,  the  transmitter  and  receiver  (TX/KX)  are  co-located. 

An  image  transmitter  TX'  at  twice  the  range  r^  from  the  receiver  leads  to  a  power 
level  equal  to  that  received  from  the  plate 


I’ 


c;.„ 


AF 


HI1 


4  IT  (2  r  ) 
P 


<n 


lil 


where 


P^,  transmitted  power, 

G^,  transmit  antenna  gain,  and 
Ap,  effective  receiving  aperture. 


TX/RX  TX/RX 


Figure  7.  Image  Source 
Calibration  Method 


TX' 


It  is  assumed  here  that  the  reflecting  plane  is  perpendicular  to  the  beam  direction 
as  indicated  in  Figure  7.  Second,  the  return  from  snow  is  determined.  With  the 
antenna  beams  directed  at  an  angle  9  (grazing  angle)  against  the  snow  surface,  the 
area  A  is  illuminated  bv  the  transmitting  antenna.  For  smell  dg.  dg  the  area  A 
is  of  elliptical  shape  and  of  size 


4  sin  9 


(2) 


with 


d...  du  3-dB  fc>  and  11-plane  transmit  antenna  beamwidths  and 
L  ri 

r  range  to  snow  surface. 

The  received  power  from  the  snow  becomes 


PR  = 


PT  GT  A  AK 
(4  Tt  rV 


(3) 


where  <7°  is  the  backscatter  cross-section  of  snow  per  unit  area. 


Note  that  Eq.  (1)  describes  the  one-way  path  loss  of  propagation  in  free  space, 
while  Eq.  (3)  describes  the  incoherent  radar  return  from  a  target  of  size  A.  Eq.  (3) 
is  simplified  from  its  original  integral  form  by  assuming  that  G^.,  o° ,  and  r  are 
constant  over  the  region  A  of  integration.  At  extremely  shallow  angles  8  this  may 
not  be  a  sufficiently  good  approximation.  The  backscatter  coefficient  0“  is  obtained 
by  combining  Eqs.  (1)  through  (3) 


4  sin  0 


<J°  = 


HP 


rr(i) 


(4) 


Only  the  received  power  ratio,  the  range  ratio,  the  grazing  angle  on  the  snow,  and 
the  transmitting  antenna  beamwidths  enter  into  the  equation.  The  receiving  antenna 
parameters  do  not  affect  the  determination  of  o°.  The  calibration  metal  plate  was 
attached  to  the  underside  of  the  sample  holder  during  the  reference  measurement, 
in  order  to  avoid  disturbing  the  snow  sample.  The  minor  range  difference  between 
metal  and  snow  surfaces  was  accounted  for  in  the  calculation.  It  is  assumed  that 
most  of  the  snow  scatter  occurs  close  to  the  surface. 

Bistatic  scatter  coefficients  are  determined  in  a  similar  wav.  Hanges  r  to  the 
transmitter  and  receiver  are  identical  by  design  and  the  same  as  in  the  backscatter 
case.  Two  beam  angles  have  to  be  considered  now:  (9  .j.  between  transmitted  beam 
and  snow  surface,  and  8  between  received  beam  and  snow  surface,  lie-examining 
Eqs.  (1)  through  (4)  one  finds  that  these  still  hold  in  the  bistatic  case  if  8  ^ .  is  sub¬ 
stituted  for  9  .  The  beamwidths  d,.  and  d,,  are  transmitter  beamwidths  as  before. 

E  II 

The  bistatic  scatter  coefficient  then  becomes 

a.  .  -  PR4°in9T  .  /  Jl\2  •  (5) 

9  T  9  R  PKP  fiE  ,!II  V  ‘  p  / 


The  attenuation  coefficient  of  snow  may  be  defined  from  a  loss  measurement 
when  the  signal  passes  through  a  layer  of  snow.  With  the  equipment  configured  as 
in  Figure  3,  power  is  received  through  the  snow.  Bather  than  making  a  refer¬ 
ence  measurement  by  removing  the  snow  slab,  the  backscatter  reference  measure¬ 
ment  is  used  as  in  Figure  1  with  the  metal  plate  mounted  under  the  basket.  For 
slab  thickness  s  we  define  the  attenuation  coefficient  as 

L  =-10  log  |  rr-~  |  '  s  *  in  dB/m  .  (6) 

V  RP/ 

This  procedure  avoids  the  inconvenience  of  removing  the  snow  slab  and  the  danger 
of  breaking  it  in  the  process. 
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In  any  of  tile  measurements,  only  one  microwave  parameter,  the  ratio 

recw*ve^  power  from  the  snow  to  received  power  from  a  reference 
target,  is  required.  Transmitted  power  levels  are  of  no  importance  and  were  not 
determined  accurately.  In  order  to  establish  the  ratio  l’^/l’^j,  in  an  individual 
measurement,  the  data  recorded  by  the  X \  -  recorder  were  calibrated  with  a  pre¬ 
cision  IT  calibrator  (III’  355  B)  inserted  after  the  first  IT  amplifier.  The  accuracy 
of  this  attenuator  was  checked  over  a  ,'O-dB  range  by  comparison  with  a  calibrated 
signal  from  a  signal  generator  (HP  tiOUA).  The  deviation  was  less  than  ±  0.  2  dll 
over  the  whole  range.  Assuming  that  the  first  conversion  in  the  receiver  is  in¬ 
herently  linear,  this  would  also  be  the  accuracy  expected  for  the  received  nun-wave 
signal.  To  check  this  assumption,  measurements  were  made  directlv  at  35  and 
1)8  GHz.  where  50-dB  variable  precision  attenuators  were  available.  The  two  trans¬ 
mitters  were  placed  2  m  from  the  receiver  with  the  variable  attenuators  connected 
between  the  sources  and  the  antennas.  Calibration  points  were  obtained  by  compar¬ 
ing  the  receiver  oul put  level  after:  (a)  inserting  a  specific  attenuation  into  Hie 
mm-wnve  attenuator,  and  (b)  replacing  it  with  the  same  attenuation  at  IT  frequency. 
Curves  are  plotted  in  Figure  8  of  the  recorder's  5  -excursion  as  a  function  of  the 
ItT  and  IT  attenuator  settings.  Note  that  two  different  curves  are  shown  for  the 
IT  attenuator  at  35  and  1)8  GHz.  This  can  be  explained  in  the  following  way.  Differ¬ 
ences  in  transmitted  power,  antenna  gain,  and  mixer  conversion  loss  lead  to  differ¬ 
ent  output  levels  at  the  XV- recorder  for  0-dB  RT  and  IT  attenuation.  During 
calibration  measurements,  the  recorder  output  was  set  to  full-scale  at  either 
frequency  by  changing  the  recorder  dc  gain.  This,  in  combination  with  the  near- 
logarithmic  receiver  response,  leads  to  different  IT  attenuator  curves  in  Figure  8. 
Consequently,  comparisons  are  possible  only  between  curves  taken  at  tile  same 
frequency. 

/It  35  GHz,  any  particular  receiver  output  amplitude  corresponds  to  IT  and  HT 
attenuator  settings  different  by  no  more  than  1  dB.  When  set  to  the  same  value, 
the  RF  attenuator  produces  an  output  lower  than  the  IT  attenuator  with  the  exception 
of  the  very  high  attenuations,  where  some  leakage  around  the  RT  attenuator  was 
noticed.  In  contrast,  at  98  GHz  the  RT  attenuator  curve  is  always  above  the  11- 
attenuator  curve,  with  the  difference  at  any  particular  receiver  output  amplitude 
not  exceeding  1  dB  in  the  attenuation  range  from  0  to  30  dll  and  not  exceeding  3  dB 
for  the  remainder  of  the  range.  Because  of  deviations  in  opposite  directions  at 
35  and  98  GHz,  it  is  assumed  that  these  are  basically  due  to  RT  attenuator  calibra¬ 
tion  error.  The  absolute  error  is  in  line  with  the  2-percent  accuracy  quoted  for 
precision  RT  attenuators.  The  increasing  deviation  at  98  GHz  in  the  range  between 
30  and  50  dB  is  attributed  to  leakage.  Leakage  does  not  affect  the  snow  measure¬ 
ments,  since  no  RT  calibration  is  performed  there.  System  calibration  at  all  fre¬ 
quencies.  including  140  GHz,  is  normally  achieved  by  variable  IT  attenuation. 
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IK  attenuation  is  known  with  better  accuracy,  and  the  35-CHz  curves  in  Figure  8 
suggest  good  KK  to  IK  conversion  linearity  over  a  50-dll  range. 


Figure  8.  Comparison  of  KK  and  IK  Calibration 


2.3  Signal  Characteristics 

Krozen  metamorphic  snow  that  can  be  picked  up  for  measurement  consists  of 
several  discernible  structures.  Snow/ice  pellets  in  various  stages  of  transition 
closely  packed  and  bonded  together  represent  the  predominant  structure.  Depending 
on  the  melt-freeze  cycle  and  on  rainfalls,  horizontal  layers  appear  of  varying 
thickness.  Sometimes  the  layers  are  of  solid  ice.  The  dielectric  properties  of  ice 
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at  0°C  have  been  given  in  a  survey  paper  by  Ray.  * 1  At  35,  98,  and  140  GHz,  the 

-4 

real  part  of  the  refractive  index  is  n^  =  1.  78.  The  imaginary  part  is  n.  =  9  X  10 
at  35  GHz,  n.  =  5  X  10  ^  at  98  GHz,  and  l  =  3X  10  ^  at  140  GHz.  Frozen  snow  is 
an  ice-air  mixture  whose  dielectric  properties  can  be  computed  from  the  fractional 

volumes  occupied  by  the  constituents.  In  an  earlier  report,  we  determined  for 

2  2 
freshly  fallen  snow  n^  =  1.  12  and  for  late  metamorphic  snow  np  =  1.  92.  l-'rom  the 

low  imaginary  component  of  the  refractive  index,  it  follows  that  any  appreciable 
attenuation  in  the  medium  is  predominantly  due  to  scattering.  Based  on  a  multitude 
of  closely  packed  low-loss  spherical  scatterers,  we  expect  to  receive  a  Rayleigh- 
distributed  signal  of  random  amplitude,  as  the  snow  slab  rotates  relative  to  trans¬ 
mitter  and  receiver  directions.  To  the  extent  that  the  ice  layers  are  capable  of 
specular  reflection,  we  expect  at  least  a  fraction  of  the  received  signal  to  show  a 
behavior  similar  to  that  signal  reflected  from  the  metal  calibration  plate. 

A  theoretical  analysis  for  densely  packed  low-loss  Alie  scatterers  does  not 
exist  at  this  time.  Lacking  a  rigorous  theory,  we  present  somewhat  qualitative 
models  in  order  to  explain  experimentally  observed  signal  amplitudes  and  fluctua¬ 
tion  rates  as  a  function  of  grazing  angle  9  .  Consider  in  Figure  9(a)  a  bistatic 
system  where  transmitter  TX  and  receiver  RX  with  angular  spacing  9  and  range 
f  are  connected  through  scatterers  Sj  and  S2>  Both  scatterers  are  within  the  snow 
slab,  Sj  at  the  center  of  rotation  and  S,,  at  a  distance  d  from  the  center,  where  d 
is  chosen  to  be  the  distance  to  the  edge  of  the  3-dB  transmit  antenna  cone.  At  this 
extreme  point  illuminated  by  the  antenna  and  for  some  angular  change  dff  of  the 
slab.  S2  travels  the  greatest  distance  of  all  scatterers  within  the  beam  and  hence 
changes  pathlength  f  j+  C2  the  most.  The  change  in  pathlength  for  a  signal  scattered 
from  S2  relative  to  the  signal  from  Sj  over  constant  path  2f  gives  rise  to  a  phase- 
shift  between  the  two  components  of  the  received  signal.  There  may  be  constructive 
or  destructive  interference.  Depending  on  the  wavelength  used  and  other  geometri¬ 
cal  parameters,  an  angle  change  d$  can  be  calculated,  required  to  change  the  phase 
angle  between  the  two  phasors  from  0  to  180  deg.  We  determine  first  the  rate  of 
range  change  as  a  function  of  grazing  angle  8  .  According  to  the  law  of  cosines,  we 
find  from  Figure  9(a) 

Cj  =  f2  +  d2  -  2  fd  cos  (9  +  7!  "  9  q)  and 

f2  =  i2  +  d2  -  2  £d  cos  8  . 


11.  Ray,  P.  S.  (1972)  Broadband  complex  refractive  indices  of  ice  and  water, 

Appl.  Optics  H(No.  8):  1836- 1844. 

12.  Lammers,  U.H.  W.  ,  and  Hayes,  D.T.  (1980)  Multipath  Propagation  Over  Snow 

at  Millimeter  Wavelengths,  In-House  Hep.  ftAl!)C-T R-80-54,  AD  A087747. 
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Figure  9.  Signal  Path  Geometries,  (a)  Bistatie  Scatter,  (b)  Reflection, 
(c)  Attenuation 


We  differentiate 

f  l  +  (  2  =  (£2  +  d2  +  2  td  cos  (0-0  Q))1/2  +  U2  +  d2  -2  fd  cos  6  )1/2  , 
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The  backscatter  case  follows  from  the  general  bistatic  case  by  setting  0  0  to  180  deg 
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The  transmission  case  (0  Q  =  0  deg,  transmitter  and  receiver  on  a  straight  line) 
yields 
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_ Cd  sin  0 _ 

(C2  +  d2  +  2  Cd  cos  0  )1/2 


(9) 


Results  calculated  from  Bqs.  (7),  (8),  and  (9)  are  plotted  in  Figure  10  and  labeled 
"Bistatie  Scatter",  "Backscatter",  and  "Attenuation  (Single  Scatter)”,  respectively. 
For  grazing  angles  0  between  90  and  15  deg,  the  change  in  pathlength  in  mm/deg 
of  slab  rotation  for  the  signal  scattered  at  S2  can  be  read  from  Figure  10.  For  each 
of  the  wavelengths  (8.6  mm,  3.0  mm,  and  2.  1  mm)  the  phase  change  for  the  signal 
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scattered  at  is  calculated  from  this  value.  A  rough  measure  is  obtained  of  how 
many  signal  fluctuations  can  be  expected  over  the  snow  slab  angular  range  of  motion. 
Note  that  only  phasing  between  scatterers  at  extreme  and  center  points  has  been 
considered. 


90  60  30  o 

GRAZING  ANGLE  S(OEG) 

Figure  10.  Diagram  to  Compute  Signal  Phasing  Rates 


In  the  backscatter  case,  pathlength  changes  of  4.  2  mm  at  0  =90  deg  and 
1.  2  mm  at  0  -  15  deg  result  from  a  1-deg  slab  rotation.  This  translates  into  a 
176-deg  phase  change  at  8.6  mm  wavelength  and  a  720-deg  phase  change  at  2.  1  mm 
wavelength  (both  at  grazing  angle  0  =  90  deg).  Conversely,  we  calculate  that  in 
order  to  go  from  constructive  to  destructive  interference  at  this  slab  orientation 
(phase  change  of  180  deg)  the  basket  must  rotate  1  deg  at  8.6  mm  wavelength  and 
0.  25  deg  at  2.  1  mm  wavelength.  At  a  grazing  angle  0  =  15  deg.  phase  changes  are 
about  3.  5  times  less  per  1-deg  basket  rotation.  Only  one  curve  of  bistatic  scatter 
is  shown  in  Figure  10.  For  0  q  =  90  deg,  the  rate  of  range  change  as  a  function  of 
slab  angle  peaks  at  6  =45  deg  at  a  level  (3  mm /deg)  comparable  to  that  found  for 
backscatter. 

Experimental  attenuation  measurements  on  snow  slabs  frequently  show  substan¬ 
tial  and  rapid  phasing  as  a  function  of  the  grazing  angle.  This  is  not  surprising  since 
scattering  from  low-loss  dielectric  spheres  should  produce  substantial  forward 
scatter.  The  curve  labeled  "Attenuation  (Single  Scatter)"  gives  phase  differences 
to  be  expected  between  signals  propagating  along  the  straight  transmitter/ receiver 
path  (0=0  deg)  and  that  path  including  scatterer  at  radius  d.  When  the  basket  is 
perpendicular  to  the  straight  line  connecting  the  transmitter  and  the  receiver,  the 
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rate  of  change  is  zero.  It  reaches  a  maximum  of  0.25  mm,  deg  at  around 
9  =  45  deg.  This  value  is  much  lower  than  the  one  for  backscatter  and  does  not 
account  for  the  rate  of  phase  fluctuations  actually  observed  experimentally.  In 
order  to  account  for  effects  observed,  three  other  interference  phenomena  involving 
the  snow  slab  have  been  investigated.  In  the  first  case  (depicted  in  .Figure  0(b))  it 
is  assumed  that  the  homogeneous  plane-parallel  dielectric  snow  slab  (nf  =  1.  99) 
refracts  and  reflects  the  microwave  signal.  Interference  is  set  up  between  the 
direct  signal  (refraction  only)  and  the  reflected  signal  (refraction  and  reflection). 

As  before,  a  calculation  of  phase  relationship  is  made.  No  consideration  is  given 
to  the  magnitude  of  the  reflection  coefficient  or  the  losses  inside  the  medium. 
According  to  Snell's  law  we  write 

sin  (90  -  9  )  =  n  sin  9  or  6  =  . 

r  r  r  n 

r 

The  path  difference,  which  we  are  interested  in,  is 


r  =  2  s,  -  s, 


Quantity  r  can  be  expressed  in  terms  of  the  slab  thickness  s.  refractive  index  n^, 
and  angle  9  shown  in  Figure  9(b) 
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Differentiation  of  Eq.  (10)  with  respect  to  9  yields 


dr 
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r  r 


2  s  ( 1  -  cos“  din)  sin  0  cos  0 

- 27; - 277  Xjn - 

n  ( 1  -  cos  9  /n  ) 
r  r 


(11) 


This  equation  is  also  plotted  in  F  igure  10  and  labeled  "Keflection ",  Typicall  ,  it 
produces  minimum  phase  change  as  the  basket  moves  through  th<  region  where  the 
direct  signal  penetrates  the  slab  perpendicularly.  We  see,  however,  that  for 
angles  8  near  40  deg,  phasing  is  as  rapid  as  for  the  backscatter  case.  The  maxi¬ 
mum  rate  of  range  change  in  this  area  is  2.6  mm /deg. 

Another  mechanism  conceivable  to  introduce  rapid  phase  fluctuations  into  the 
attenuation  measurement  is  termed  "Attenuation  (Multiple  Scatter)"’  in  Figure  10. 

We  assume  the  existence  of  at  least  dual  scattering  as  sketched  in  Figure  9(c).  Here 
the  direct  transmitter- receiver  path  establishes  one  signal  phase  and  the  path 
through  Sj  and  Sg  the  second  one.  If  this  propagation  mode  is  a  valid  one,  inspection 
of  Eq.  (9)  shows  the  appropriate  expression 
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Eq.  (12)  can  account  for  rapid  phasing  in  an  attenuation  measurement,  actually  at 
half  the  rate  observed  for  backscatter.  This  is  apparent  from  an  inspection  of 
Figure  10. 

A  last  interference  mechanism  affecting  the  attenuation  measurement  is  poss¬ 
ible  by  diffraction.  At  low  grazing  angles  6  ,  radiation  diffracted  by  the  edge 
of  the  basket  toward  the  receiver  may  be  of  comparable  strength  to  the  radiation 
transmitted  directly  through  the  snow  slab.  This  is  especially  possible  if  the  direct 
signal  has  been  heavily  attenuated  by  the  snow.  In  terms  of  the  rate  range  change, 
this  mechanism  is  equal  to  single  scatter  as  described  in  Eq.  (9).  The  only  differ¬ 
ence  lies  in  the  distance  d  from  the  center  of  rotation.  If  d  =  0.  .'18  m  is  used  for 
the  edge  of  the  basket  rather  than  d  =  0.  12  m  for  the  edge  of  the  3-dB  beam,  a 
curve  in  Figure  10  results  quite  similar  to  that  labeled  "Reflection'1.  There  is  a 
substantial  increase  in  the  rate  of  range  change  over  that  for  "Attenuation  (Single 
Scatter)".  This  more  rapid  phasing  should  be  noticeable  only  when  the  basket  edge 
rotates  enough  to  approach  the  edge  of  the  transmit  beam.  Experimental  results 
in  Section  4  will  be  compared  against  those  shown  in  Figure  10. 

The  change  in  received  power  as  a  function  of  grazing  angle  is  another  param- 

1 3 

eter  aiding  in  the  interpretation  of  millimeter- wave/ snow  interaction.  Clapp,  in 
an  early  report  on  terrain  reflection  as  seen  by  airborne  radar,  discusses  three 
simple  models:  (a)  diffuse  reflection  by  the  rough  surface  of  a  lossless  medium, 

(b)  omnidirectional  direct  reflection  from  a  single  layer  of  spherical  particles  with 
random  mutual  spacing,  located  approximately  at  equal  height  above  a  plane  re¬ 
flecting  surface,  and  (c)  omnidirectional  reflection  from  the  spherical  particles  of 
multiple  stacked  layers  as  in  (b).  such  that  multiple  scattering  occurs.  The  depth 
of  layers  is  large  enough  to  render  background  contributions  to  the  backscattered 
signal  insignificant.  The  absorption  factor  (1-6)  as  well  as  the  spatial  arrange¬ 
ment  of  spheres  and  the  grazing  angle  relative  to  the  "surface"  of  the  medium 
determine  the  relative  importance  of  first  and  higher  order  reflections. 

Models  (a)  and  (c)  have  the  potential  to  represent  granular  snow  quite  well. 

Clapp  assumes  in  (a)  an  angular  distribution  of  scattered  power  according  to 
Lambert's  law.  For  unidirectional  illumination  as  in  the  radar  case,  this  can  be 
achieved  by  a  dissipationless  dielectric.  Due  to  multiple  reflections  and  refractions, 
its  reradiation  characteristic  is  independent  of  the  direction  of  the  incident 


13.  Clapp,  R.  E.  (1946)  A  Theoretical  and  Experimental  Study  of  Radar  Ground 
Return,  Rep.  1024  Rad.  Lab.  MIT,  Cambridge,- Massachusetts. 
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radiation.  Keradiated  power  density  changes  with  the  sine  of  grazing  angle  9  . 

Since  the  incident  power  per  unit  surface  area  is  also  proportional  to  sin  8  ,  back- 
scatter  from  a  rough  surface  obeying  I.ambert's  law  is  proportional  to  sin^f?  . 

Model  (c)  is  based  on  the  assumption  of  a  mean  free  pathlength  between  spheres 
much  smaller  than  the  thickness  of  the  medium.  This  is  the  case  for  typical  snow 
slabs.  Another  assumption  is  that  of  uniform  reradiation  from  the  spheres.  This 
may  or  may  not  be  true  for  granular  ice  particles  at  millimeter  wavelengths.  For 
lack  of  a  better  model  it  will  have  to  be  accepted  here.  If  the  spheres  are  lossy 
(6  -*  0),  single  reflections  predominate,  leading  to  an  angular  dependence  of  re¬ 
flected  power  proportional  to  sin  8  .  With  6  -*  1,  second,  third,  and  higher  order 
reflections  contribute  significantly  to  the  backscattered  power.  The  0  -dependence 
becomes  more  complex.  Equations  based  on  up  to  third  order  reflections  are  in¬ 
cluded  in  Clapp's  report.  Kor  the  purpose  of  comparison  with  our  experimental 
data  we  have  replotted  in  Figure  11  his  model  (a)  rough  surface  scatter  and  model 
<c)  volume  scatter  including  third  order  reflections,  both  for  a  lossless  medium. 

In  addition,  the  specular  reflection  model  for  backscatter  is  shown  in  the  same 
figure.  It  is  simply  determined  by  the  product  of  transmit  and  receive  beam  pat¬ 
terns.  Since  actual  3-dB  beamwidths  range  from  6.  8  to  13.  9  deg,  the  curve  in 
Figure  1 1  is  for  an  average  9  deg.  All  three  curves  have  been  normalized  to  0  dB 

at  perpendicular  incidence  ( 8  =  90  deg).  The  curve  labeled  "Rough  Surface  Scatter" 
2 

follows  sin  6  exactly,  while  the  curve  labeled  "Volume  Scatter"  drops  off  more 
strongly  with  decreasing  8  than  sin  8  would.  This  is  the  effect  of  multiple  scatter¬ 
ing.  Note  that  the  increase  in  scattered  signal  intensity  in  the  perpendicular  direc¬ 
tion  ( 9  =  90  deg)  due  to  multiple  reflections  has  been  eliminated  in  the  normaliza¬ 
tion  process. 


Figure  11.  Angular 
Dependence  of  Backscattered 
Power 
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2.4  Measurement  Procedure:  and  1. oca t ions 


One  area  of  concern  in  the  calibration  of  the  sample  measurements  is  the  re¬ 
turn  from  the  waxed  string  or  nylon  fishing  line  holding  the  snow,  and  from  the 
wooden  frame  of  the  sample  holder.  A  series  of  measurements  were  made  with 
the  sample  holder  empty  but  fully  strung  with  erosswoven  5  >  5-cm  netting,  and 
also  with  the  sample  holder  empty  and  some  or  all  of  the  strings  removed.  The 
majority  of  data  in  this  report  were  obtained  with  all  strings  mounted  along  the 
direction  perpendicular  to  the  basket  axis  of  rotation  and  two  each  at  the  end  of  the 
basket  crosswise  for  better  support  of  the  slab  (Figure  1).  Measurements  on  the 
empty  holder  were  made  for  horizontal,  vertical,  or  crossed  polarization  at  the 
three  frequencies.  Returns  from  the  netting  relative  to  that  from  the  horizontally 
or  vertically  polarized  return  from  the  metal  plate  are  listed  in  Table  2. 


Table  2.  Uaekscatter  From  Sample  Holder 


Polarization 

!• 

35 

requenev  (GHz) 
98 

140 

HH 

-50  dll 

-30  dll 

-38  dB 

VV 

-50  dB 

-2G  dB 

-29  dll 

HY 

-50  dll 

-50  dll 

-48  dll 

Here,  H  stands  for  horizontal  and  V  stands  for  vertical  polarization,  with  the  first 
letter  denoting  the  transmitted  and  the  second  letter  denoting  the  received  polariza¬ 
tion.  In  none  of  the  snow  slab  measurements  was  the  return  less  than  -20  dll.  The 
contribution  of  the  netting  is  therefore  considered  insignificant. 

Snow  measurements  during  the  winter  of  1977/78  were  closely  linked  to  the 

evolution  of  the  equipment  from  the  sample  configuration  to  the  in-situ  configura- 

14 

tion.  Measurements  during  the  early  part  of  the  winter  were  conducted  at  the 
Prospect  Hill  Field  Station  at  an  elevation  of  145  m  above  sea  level.  1'he  equip¬ 
ment  was  positioned  on  a  loading  platform  with  a  major  part  of  it  inside  a  heated 
building  and  only  the  boom  extending  outside.  A  plastic  window  closed  the  opening 
and  prevented  a  temperature  change  around  the  snow  slab  due  to  heat  generated 
inside  the  building.  The  loading  platform  afforded  an  additional  1-m  separation  of 
the  test  path  from  the  ground,  further  reducing  the  possibility  of  unwanted  reflection. 


14.  Ilayes,  D.T.,  Lammers,  U.il.W.,  and  Marr,  R.A.  ( 1 98 1 1  Scattering  F  rom 
Snow  Background  at  35,  98,  and  140  Clilz,  KALX'-TR  (to  be  published). 


Snow  slabs  of  proper  size  were  cut  from  undisturbed  areas  in  the  vicinity.  In 
most  cases  slabs  were  placed  in  the  sample  holder  and  left  there  until  completion 
of  the  measurements.  A  complete  set  of  measurements  on  a  particular  sample  con¬ 
sisted  of  HH,  VV,  and  HV  scans  at  all  frequencies,  including  a  calibration  against 
the  metal  plate  before  and  after  each  scan.  Sample  measurements  were  begun 
16  January  1978  and  concluded  2  March  197  8.  During  this  time,  data  from  a  total 
of  230  runs  were  collected. 


3.  SNOW  CHARACTERIZATION 

Snow  deposited  on  the  ground  undergoes  substantial  changes  while  aging  under 
various  weather  conditions.  Changes  in  the  physical  properties  of  snow  cover  have 
a  profound  effect  on  its  interaction  with  millimeter  waves.  An  effort  was  therefore 
made  during  the  experimental  program  to  determine  in  detail  the  relevant  physical 
properties  of  snow  at  or  near  the  location  where  samples  were  removed. 

3.1  Snow  Development 

In  one  aspect,  the  winter  of  1977/78  was  typical  for  the  Boston  area.  The 
major  snowstorms  took  place  during  the  months  of  January  and  February.  The 
winter  was  atypical  in  the  cumulative  amount  of  snow  deposited  and  the  large 
amounts  delivered  by  several  storms.  A  total  of  215  cm  of  snow  fell  in  the  Boston 
area.  This  was  more  than  double  the  av  rage  snowfall  (100  cm)  expected  in  a  single 
season.  On  two  occasions,  record  snowstorms  occurred.  The  storm  of  20  January 
deposited  56  cm  of  snow.  This  was  a  record  for  a  24 -hr  period.  On  6  February 
the  "biggest  snowstorm  ever"  began.  By  the  time  this  storm  ended  on  7  February, 
76  cm  of  snow  had  been  deposited.  High  winds  associated  with  the  storm,  coming 
at  the  time  of  full  moon,  resulted  in  severe  coastal  flooding  and  destruction.  At 
the  conclusion  of  this  storm  an  accumulated  total  of  180  cm  had  been  deposited 
since  the  first  snowstorm  on  26  November  197  7.  Deep  snow  in  the  area  generally 
does  not  last  for  a  long  time.  Boston  s  location  on  the  New  England  coast  ensures 
a  variability  of  weather  (cold  to  warm,  snow  to  rain)  to  prevent  any  buildup  of  snow 
over  a  full  season.  When  measured  a  week  after  the  6  February  storm,  the  snow 
depth  at  the  Prospect  Hill  Field  Station  was  86  cm.  Of  this,  61  cm  came  from  the 
last  storm.  The  remainder  had  been  deposited  during  two  previous  storms. 

Whether  the  Boston  area  is  to  be  hit  by  snow,  rain,  or  a  mixture  of  the  two  de¬ 
pends  upon  the  track  of  the  low  pressure  area  which  approaches  from  the  south.  If 
the  low  tracks  to  the  west  of  the  city,  rain  is  brought  on  southwesterly  winds.  If 
there  is  a  large  mass  of  cold  air  hanging  over  Boston  as  the  low  approaches  from 
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from  the  south,  the  precipitation  begins  as  snow.  As  the  warmer  air  circulating 
around  the  low  displaces  this  colder  air,  the  precipitation  changes  to  sleet  or 
freezing  rain  and  eventually  goes  over  to  rain.  On  the  other  hand,  if  the  low  passes 
to  the  east  of  the  city,  precipitation  in  the  form  of  snow  is  the  result.  In  November, 
the  most  probable  track  of  a  storm  is  to  the  west  of  Boston.  As  the  winter  pro¬ 
gresses,  the  most  probable  path  moves  eastward.  In  January,  it  lies  to  the  east 
of  Boston.  Thus,  November  and  December  storms  are  more  likely  to  produce  rain 
or  start  off  as  snow  and  end  as  rain.  January  and  February  storms  are  most  like¬ 
ly  to  bring  only  snow.  Precipitation  in  the  form  of  rain  is  not  an  improbable  event, 
even  in  the  later  winter  months.  However,  there  is  usually  enough  snow  on  the 
ground  at  this  time  so  that  not  all  is  washed  away  if  rain  does  occur. 

It  should  be  clear  from  the  above  discussion  that  winter  weather  in  the  Boston 
area  is  most  variable  and  that  precipitation  in  the  form  of  rain  is  a  possibility 
throughout  the  winter.  Thus,  partial  melting  and  subsequent  refreezing  due  to 
temperature  changes  from  day  to  night  and  drenching  by  rain  play  an  important  role 
in  modifying  the  structure  of  the  snow  already  deposited  on  the  ground. 

Newly  fallen  snow  in  crystalline  form  will  begin  to  undergo  changes  almost 
immediately.  If  the  vertical  temperature  gradient  through  the  snowpack  is  small, 
water  molecules  tend  to  sublime  from  the  sharper  edges  of  the  snow  crystal  where 
the  vapor  pressure  is  high  and  condense  in  sections  of  the  crystal  where  the  radius 
of  curvature  is  low  (low  vapor  pressure).  This  breaks  up  the  crystalline  struc¬ 
ture  of  the  snow.  The  snow  crystals  slowly  metamorphose  into  small  round  ice 
grains.  The  ice  grains  continue  to  undergo  changes.  Vapor  sublimes  from  the 
smaller  ones  and  condenses  on  the  larger  ones.  This  entire  process  is  called  equi- 
temperature  metamorphism  (ET).  The  rapidity  of  the  change  is  proportional  to  the 
temperature  of  the  snowpack.  It  is  most  rapid  near  0°C  and  almost  ceases  at 

-40°C.  When  equilibrium  is  reached,  the  medium  consists  of  grains  of  nearly 
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identical  size,  and  the  density  of  the  snow  is  roughly  0.6  g/cm  . 

If  a  strong  temperature  gradient  exists  in  the  snow,  the  water  vapor  tends  to 
migrate  from  the  warmer  region  (near  the  ground)  toward  the  colder  region,  where 
it  condenses.  This  is  temperature  gradient  (TG)  metamorphism.  It  is  rarely  seen 
in  New  England. 

The  dominant  metamorphic  process  affecting  snow  in  the  Boston  area  is  the 
so-called  melt-freeze  (MF)  metamorphism.  Typically,  following  a  midwinter  storm, 
high  pressure  settles  over  the  area  bringing  in  cold  air.  During  this  time,  ET 
metamorphism  is  dominant.  Eventually  the  high  passes  over  the  region,  bringing 
southerly  winds.  The  air  temperature  rises  above  freezing  during  the  day  or  at 

15.  Perla,  R.  I.  ,  and  Martinelli,  M.  ,  Jr.  (1976)  Avalanche  Handbook,  Agriculture 
Handbook  489,  U.  S.  Dept,  of  Agr.  ,  Forest  Service. 
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least  high  enough  that  solar  radiation  causes  appreciable  melting  to  take  place  in 
the  upper  layer  of  the  snow.  The  high  may  be  followed  by  a  low  bringing  in  rain. 

In  any  case,  free  water  appears  in  the  upper  snow  layer.  It  filters  down  through 
the  snow,  forms  a  liquid  layer  around  a  snow  grain,  and  possibly  connects  two 
grains  by  a  liquid  meniscus.  At  night  the  water  refreezes.  As  a  net  result  the  ice 
grains  grow  in  size  (up  to  2  or  3  mm)  and  the  bearing  strength  of  the  snow  layer 
increases.  This  process  can  take  place  most  rapidly.  Late  in  the  season,  just 
one  warm  day  with  strong  solar  heating  followed  by  a  cold  night  is  enough  to  destroy 
all  vestiges  of  the  virgin  snow  and  produce  a  hard  layer  of  ice  grains. 

3.2  Snow  Measurement 

To  provide  the  necessary  supporting  data  for  analysis  of  the  electromagnetic 
scattering  and  attenuation,  selected  characteristics  of  the  snow  cover  were  re¬ 
corded  concurrently  with  the  electromagnetic  measurements.  In  addition,  the  de¬ 
velopment  of  the  snow  cover  during  the  winter  snow  season  was  recorded  in  a  daily 
weather  diary.  The  following  snow  properties  were  measured  on  a  daily  basis: 

(1)  Depth, 

(2)  Density, 

(3)  Hardness, 

(4)  Temperature, 

(5)  Stratigraphy, 

(6)  Microstructure, 

(7)  Surface  characteristics. 

The  majority  of  these  properties  were  determined  using  equipment  supplied  in  the 
U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL)  snow 
observation  kit,  manufactured  by  Geotest  Instrument  Corporation,  Wheeling,  Illinois. 
With  the  exception  of  temperature,  all  properties  were  assumed  to  remain  un¬ 
changed  during  a  daily  test  period.  The  snow  was  characterized  each  day  before 
the  electromagnetic  measurements  were  begun.  The  snow  properties  were  re¬ 
corded,  and  a  surface  sample  (4  to  14  cm  thick)  was  cut  to  be  used  in  the  scatterom- 
eter.  A  trench  running  into  an  area  of  unperturbed  snow  was  established.  Each 
day  the  trench  was  extended  exposing  virgin  snow. 

The  following  measurement  techniques  pertain  to  individual  snow  parameters: 

( 1 )  Snow  Depth 

The  total  snow  depth  in  the  trench  was  recorded  daily. 

(2)  Snow  Density 

Snow-sampling  tubes  provided  in  the  observation  kit  were  6  cm  in  diameter. 

Their  volume  was  500  cm3.  At  least  three  horizontal  samples  were  ob¬ 
tained  daily  from  each  snow  layer. 
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(3)  Hardness 


The  surface  hardness  was  measured  using  a  Canadian  hardness  gauge 
supplied  with  the  CRREL  kit.  This  gauge  consists  of  a  push-type  spring 
balance  with  provision  for  mounting  disks  of  various  sizes  on  the  end 
of  the  push  rod.  in  operation,  the  disk  is  pressed  against  the  snow  sur¬ 
face  and  the  reading  on  the  push  rod  scale  is  noted  when  the  surface 

2 

ruptures.  The  gauge  is  calibrated  in  g/cm  . 

(4)  Temperature 

Himetallic  thermometers  (range  -50  to  20°C),  supplied  in  the  CRREL  kit, 
were  used  to  obtain  a  temperature  profile  of  the  snow  cover.  At  least 
one  reading  per  layer  was  made.  For  layers  thicker  than  10  cm,  the  snow- 
temperature  was  measured  at  about  5-cm  intervals.  The  air  temperature 
was  recorded  at  the  same  time. 

(5)  Stratigraphy 

The  thickness  of  all  snow  layers  was  measured.  The  presence  of  any  ice 
layers  and  their  thickness  were  noted. 

(6)  Microstructure 

A  sample  of  snow  in  each  layer  was  examined  to  determine  grain  shape 
and  size  distribution.  Photomacrographs  of  the  snow  grains  were  made 
with  a  millimeter  grid  in  the  background  in  order  to  record  the  grain 
shape  and  size  distribution. 

(7)  Surface  Characteristics 

The  surface  properties  were  described  from  visual  observation.  Also, 
photographs  of  the  surface  were  made.  Data  on  mean  surface  height 
variation  and  surface  features  were  recorded.  Types  of  surface  features 
noted  were:  ice  layer  at  surface,  roughness  caused  by  windblown  snow, 
or  snow  melt-freeze  cycle. 

Snow  data  are  presented  in  the  appendix.  Data  are  given  for  every  day  on  which 
electromagnetic  tests  were  conducted.  Reference  will  be  made  to  these  data  when 
discussing  experimental  results  in  Section  4. 


4.  EXPERIMENTAL  RESULTS 

No  free  water  was  detectable  in  any  of  the  snow  slabs  used  for  measurement. 
The  air  temperature  never  rose  above  (PC  during  tests.  The  tests  were  always 
carried  out  in  the  shade  to  eliminate  the  possibility  of  melting  caused  by  direct  solar 
radiation.  Signal  strength  was  recorded  by  pen  and  ink  plot  only,  with  received 
amplitude  along  the  Y-axis  and  sample  holder  angular  position  along  the  X-axis. 
Typical  sample  data  are  shown  in  Figures  12  through  15. 
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4. 1  Sample  Cases  of  Scatter  and  Attenuation 

It  is  possible  to  divide  the  sample  data  period  into  two  distinctive  sections. 

The  first  group  of  data  was  collected  from  16  January  through  (i  February.  During 
this  period  the  weather  was  quite  variable.  A  number  of  snowstorms  in  curred  with 
varying  amounts  of  snowfall.  Each  storm  ended  with  sleet  or  freezing  ruin,  or  the 
snowstorm  was  followed  a  few  davs  later  bv  periods  of  modi  rate  to  lie. ivy  rum.  A 
total  of  41  backscatter  measurements  were  performed  during  this  period  including 
frequencies  35,  98,  and  140  GHz  and  polarizations  \  \  and  )l\  .  l'hrre  wen-  lit' 
attenuation  measurements  at  the  three  frequencies  and  at  i  V  polarization.  Mu 
bistatic  scatter  measurements  were  made. 

Due  to  the  heavy  blizzard  on  6/7  February,  non-essential  travel  was  restricted 
statewide.  Snowclogged  roads  leading  to  the  test  site  were  not  cleared  until 
21  February,  on  which  day  measurements  were  resumed.  The  second  period  of 
sample  measurements  included  21  February  to  2  March.  The  weather  following 
the  blizzard  was  marked  by  its  uniformity.  The  days  were  sunny  for  the  most  part 
and  the  temperature  rarely  exceeded  freezing.  Only  on  three  days  did  this  occur: 
2°C  on  16  F  ebruary,  3°C  on  17  February,  and  1  ’('  on  24  February.  The  change  in 
the  deposited  snow  crystals  proceeded  by  ET  metamorphism  at  a  very  slow  rate. 

By  the  time  sample  measurements  were  resumed  on  21  February,  the  structure  of 
the  snow  had  reached  a  fairly  stable  equilibrium  condition.  During  the  course  of 
measurements  the  properties  of  the  snow  essentially  did  not  change.  The  variability 
of  the  samples  collected  on  different  days  could  be  attributed  to  the  horizontal 
inhomogeneities  in  the  snow  cover.  These  were  caused  by  drifting  and  other  factors 
such  as  variability  in  solar  heat  loading.  Some  samples  wore  collected  from  areas 
which  were  continuously  in  the  shade.  A  total  of  7  1  backscatter  and  53  bistatic 
scatter  measurements  were  performed  during  this  time.  In  both  cases  data  were 
obtained  at  35,  98,  and  140  GHz  and  at  HH,  VV,  and  1 1 V  polarization.  Attenuation 
measurements  numbered  39  at  the  same  frequencies  and  HH  and  VV  polarization. 

From  the  large  set  of  data,  three  days  were  selected  for  discussion  in  this 
section.  The  specific  days  chosen  display  both  distinctive  physical  characteristics 
of  the  snow  and  typical  effects  on  the  electromagnetic  signals.  Figures  12  and  13 
represent  unprocessed  data  from  the  earlier  period.  They  are  made  up  of  six  in¬ 
dividual  measurements  each.  Backscatter  and  attenuation  at  vertical  polarization 
are  shown  for  the  three  frequencies.  The  days  chosen,  17  January  and  28  January, 
are  near  the  beginning  and  the  end  of  the  earlier  period.  They  represent  the  vari¬ 
ability  of  snow  conditions  and  measured  data  during  this  time. 

During  this  period  two  specific  types  of  snow  cover  existed.  On  the  first  two 
days  of  testing,  the  top  5  cm  of  snow  consisted  almost  entirely  of  spherical  sleet 
particles.  Data  collected  from  a  slab  containing  only  the  sleet  layer  are  shown  in 
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Figure  12.  Subsequent  snowstorms  covered  the  sleet  layer.  Periods  of  rain 
followed  which  produced  rapid  change  of  the  new  top  snow  and  the  sleet  layer. 
Examination  of  the  layers  showed  typical  effects  of  MF  metamorphism.  The  sleet 
layer  was  perturbed  by  the  rain  to  the  point  that  it  differed  from  the  other  layers 
only  by  a  higher  density.  As  the  snow  continued  to  age.  the  thickness  of  ice  layers 
delineating  the  boundary  between  snow  layers  increased.  The  second  set  of  back- 
scatter  and  attenuation  data  (Figure  13)  was  collected  using  a  14-cm  thick  sample 
of  this  type  of  snow. 
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Figure  12.  Sample  Backscatter  and  Attenuation  Data  of 
17  January  1378,  (vertical  polarization) 

The  third  set  is  typical  of  data  collected  after  the  blizzard.  At  this  time  the 
measurement  protocol  had  been  expanded  to  include  both  horizontally  polarized  and 
bistatic  measurements. 
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Figure  13.  Sample  Backscatter  and  Attenuation  Data  of 
28  January  1978.  (vertical  polarization) 


Figure  14  shows  original  co-  and  cross-polarized  backscatter  data  collected  at  each 
frequency  on  23  February.  Figure  15  displays  bistatic  scatter  and  attenuation  data 
collected  on  the  same  day.  For  approximately  three  weeks  following  the  blizzard, 
the  weather  was  sunny  and  cold,  causing  ET  metamorphism.  The  average  snow 
properties  present  in  the  periods  from  which  these  sample  data  were  chosen  are 
compiled  in  Table  3.  The  density  (0.  26  g/cm  )  of  the  samples  collected  following 
the  blizzard  did  not  differ  much  from  that  of  the  samples  used  in  the  latter  part  of 

3 

the  January  tests  (0.  29  g/cm  ).  However,  there  was  a  marked  difference  in  the 
bearing  strength  of  the  snow  cover  between  that  produced  by  the  MF  metamorphism 
in  January  and  that  produced  by  ET  metamorphism  following  the  blizzard.  The 
rains  in  January  provided  ample  liquid  water  which  upon  freezing  increased  the 
strength  of  the  bonds  between  adjacent  snow  grains.  This  resulted  in  a  rigid  mater¬ 
ial  from  which  samples  could  be  cut  easily  and  removed  without  destroying  the 
integrity  of  the  sample. 
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Figure  14.  Sample  Backscatter  Data  of  23  February  197  8, 
(horizontal  and  cross-polarization) 


During  KT  metamorphism,  bonding  between  adjacent  snow  grains  is  produced 
primarily  by  sublimation.  As  a  result,  the  bonds  tend  to  be  weak.  This  was  the 
case  with  the  snow  after  the  blizzard.  The  top  inch  was  a  fairly  rigid  crust  pro¬ 
duced  by  solar  melting  during  the  day  and  refreezing  at  night.  Below  this  the  snow 
was  fragile.  Samples  could  be  cut  but  great  care  had  to  be  taken  so  as  not  to  break 
the  sample. 

The  data  from  each  day  will  be  considered  in  turn.  The  ordinate  of  each  graph 
is  labeled  in  dB  below  the  signal  received  from  the  metal  calibration  plate.  The 
compressed  scale  for  lower  signal  levels  is  due  to  the  receiver  AGC  characteristic. 
The  abscissa  is  marked  in  degrees  grazing  angle  of  the  transmitted  beam  on  the 
snow  surface.  The  grazing  angle  increases  from  0  to  90  deg  and  then  decreases 
to  60  deg,  even  though  the  graph  represents  a  one-directional  motion  of  the  basket, 
covering  an  angle  of  120  deg.  This  is  done  so  that  the  grazing  angle  remains  an 
acute  angle  between  the  transmitted  beam  and  the  snow  surface.  The  backscatter 
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Figure  15.  Sample  Bistatic  Scatter  and  Attenuation  Data  of 
23  February  1978,  (horizontal  polarization) 


return  at  35  GHz  from  the  sleet  sample  (6  igure  12)  is  quite  specular  in  nature.  For 
comparison  purposes  the  return  at  35  GHz  from  the  calibration  plate  is  shown  in 
Figure  16.  The  return  from  the  sleet  slab  is  about  10  dB  below  the  peak  signal 
from  the  flat  plate.  This  is  close  to  the  reflected  power  one  would  calculate  from 
a  flat  slab  of  ice.  Using  Stratton's expression  for  the  power  reflection  coefficient 
(perpendicular  incidence  on  lossless  dielectric) 

n  -  1 

K(dB)  =  20  log  ---  +  —  ,  M3) 

1  r 

one  obtains  H  =  -11  dB  for  n  =  1.  7  8,  the  value  for  ice  given  bv  Rav. 

r 

16.  Stratton,  J.A.  (1941)  Electromagnetic  Theory,  McGraw-Hill,  New  York. 
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Table  3.  Average  Parameters  of  Snow  Cover 


Time  Period 

Average 

Density 

vg/cm3) 

Average 

Grain  Size 
(mm) 

Grain 

Shape 

Comments 

16  January 

and 

17  January 

0.  53 

2 

Spherical 

Sleet 

(Firm  structure) 

27  January 
to 

6  February 

0.  29 

i 

(Many  bonded 
together  to  2, 
some  to  4  mm) 

Irregular 

Me  It -Freeze  Snow 
(Very  firm 
structure) 

21  February 
to 

2  March 

0.  26 

1 

Irregular 

Equitemp.  Snow 
(Very  brittle 
structure) 
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Figure  16.  Backs catter  From 
Flat  Metal  Plate  at  35  GHz, 
(vertical  polarization) 


The  backscattered  signal  at  slab  directions  other  than  perpendicular  exhibits 
the  scintillations  one  expects  from  a  random  collection  of  scatterers,  in  this  case 
the  sleet  particles. 

At  98  and  140  GHz  much  less  of  a  specular  response  is  seen.  The  appearance 
of  the  fluctuating  signal  is  essentially  unchanged  at  all  grazing  angles.  At  both 
frequencies,  the  peak  return  at  perpendicular  incidence  is  again  below  that  from 
the  flat  plate  by  approximately  10  dB.  Thus  it  is  roughly  equal  to  the  power  re¬ 
flected  from  an  ice  slab.  As  the  grazing  angle  decreases,  the  backscatter  signal 
decreases  only  slightly.  It  must  be  assumed  that  an  effective  backscatter  mecha¬ 
nism  exists  within  the  invididual  scattering  particle.  Incoherent  contributions  from 
many  particles  add  up  to  a  strong  sum  signal.  Backscatter  from  dielectric  spheres 
near  resonance  could  account  for  this  result. 

The  character  of  the  attenuation  data  also  differs  between  35  GHz  and  the  two 
higher  frequencies.  With  the  exception  of  grazing  angles  below  30  deg,  no  attenua¬ 
tion  was  measured  at  35  GHz.  The  transmitted  signal  is  independent  of  the  angular 
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position  of  the  basket.  At  the  two  higher  frequencies  the  transmitted  signals  are 
strongly  attenuated  and  modulated  in  a  complex  fashion.  The  modulation  of  the 
transmitted  signal  at  angles  larger  than  30  deg  appears  to  be  composed  of  a  slowly 
varying  component  of  large  amplitude  with  a  more  rapidly  varying  component  of 
lower  amplitude  superimposed  on  it.  Below  30  deg  there  is  a  change  in  the  charac¬ 
ter  of  the  modulation. 

A  number  of  changes  can  be  seen  in  the  data  (Figure  13)  collected  from  melt- 
freeze  snow.  First,  there  is  no  trace  of  a  specular  return  at  35  GHz.  Second, 
the  backscattered  signal  at  lower  grazing  angles  is  higher  than  the  return  measured 
from  the  sample  containing  sleet.  Third,  the  35-GHz  signal  is  attenuated  and 
modulated  upon  transmission  through  the  snow.  The  differences  in  the  nature  of 
the  35-GHz  data  collected  from  the  two  samples  allow  the  following  conclusion.  In 
the  case  of  the  slab  containing  sleet,  the  strong  coherent  specular  reflection  and 
the  strong  coherent  transmitted  signal  completely  dominate  any  incoherent  scatter¬ 
ing  caused  by  individual  snow  grains.  Only  for  backscatter  at  grazing  angles  such 
that  the  receiving  antenna  does  not  intercept  the  specularly  scattered  power  does 
one  see  a  fluctuating  signal  indicative  of  scattering  from  a  collection  of  random 
scatterers. 

In  contrast,  the  sample  of  melt-freeze  snow  causes  strong  incoherent  scatter¬ 
ing.  The  coherent  backward  and  forward  propagating  signals  are  decreased  to  the 
point  that  incoherently  scattered  power  dominates  and  one  sees  a  randomly  modu¬ 
lated  signal  at  all  grazing  angles.  The  backscatter  signals  at  the  higher  frequencies 
are  somewhat  lower  than  those  for  the  corresponding  sleet  data.  Aside  from  this, 
the  character  of  the  backscatter  signal  is  not  changed.  The  transmitted  signals  are 
different.  First,  they  are  more  strongly  attenuated.  This  is  to  be  expected  as  the 
slab  is  thicker  (14  cm  as  opposed  to  the  5  cm  sleet  slab).  Second,  the  modulation 
frequency  of  the  transmitted  signal  is  higher.  At  140  GHz,  the  fluctuation  frequency 
of  the  transmitted  signal  is  comparable  to  that  of  the  backscattered  signal.  This 
again  may  be  due  to  the  thicker  slab  giving  more  opportunity  for  multiple  scatter¬ 
ing.  Also,  individual  ice  particles  in  melt-freeze  snow  may  be  stronger  scatterers. 
V  isual  examination  of  the  snow  showed  clusters  of  ice  grains  forming  macrograins, 
roughly  spherical  in  shape  with  diameters  as  large  as  1  cm.  Herman  and  Battan'  1 
have  calculated  the  backscatter  cross  section  of  ice  spheres.  They  have  shown  that 
the  backscatter  cross  section  (normalized  to  the  geometrical  cross  section)  is  an 
oscillatory  function  whose  mean  value  increases  as  a  function  of  2na/A  (a  is  the 
sphere  radius).  It  reaches  a  n  aximum  of  40  at  2na/A  =  60,  before  decreasing 
asymptotically  to  0.078  as  2na/A  approaches  infinity.  Increased  snow  grain  size 
could  produce  greater  backscatter  and  attenuation. 

17.  Herman,  B.  M.  ,  and  Battan,  L.  J.  (1961)  Calculations  of  Mie  backscattering  of 
microwaves  from  ice  spheres,  Quart.  J.  of  the  Roy.  Met.  Soc.  87:223-230. 


Finally,  consider  the  unprocessed  data  collected  on  23  February.  A  specular 
backscatter  signal  is  seen  at  35  GHz  (left-hand  side  of  Figure  14),  but  it  appears 
weaker  than  that  measured  off  the  sleet  sample  on  17  January.  The  backscatter 
returns  determined  at  98  and  140  GHz  are  very  similar  to  those  from  the  sleet 
and  melt-freeze  snow  samples.  The  cross-polarized  backscatter  signal  is  shown 
on  the  right-hand  side  of  Figure  14.  The  35-GHz  cross-polarized  return  displays 
no  specular  response  and  is  always  more  than  35  dB  below  the  co-polarized  rnetal- 
plate  return.  Compared  to  the  co-polarized  signal  from  the  snow,  the  cross- 
polarized  signal  was  about  20  dB  below  at  perpendicular  orientation  and  about  10  dB 
below  at  other  angles.  The  situation  is  different  at  98  and  140  GHz.  At  98  GHz  the 
cross-polarized  return  is  never  more  than  4  dB  below  the  co-polarized  return. 

The  co-  and  cross-polarized  signals  are  somewhat  comparable  also  at  140  GHz. 

At  both  frequencies  the  appearance  of  the  fluctuating  signals  is  similar. 

The  bistatic  return  is  shown  on  the  left-hand  side  of  Figure  15.  During  the 
bistatic  measurements,  the  difference  between  the  transmit  and  the  receive  beam 
directions  was  set  to  90  deg.  All  bistatic  data  were  obtained  at  this  single  constant 
scatter  angle.  At  each  frequency  the  envelope  of  the  fluctuating  signal  has  a  maxi¬ 
mum  at  45  deg  grazing  angle  and  falls  off  to  either  side  of  this  angle.  At  45  deg, 
we  expect  the  maximum  specularly  reflected  signal  to  occur.  The  drop-off  is 
sharpest  at  35  GHz.  The  specular  appearance  at  this  frequency  is  comparable  to 
that  in  the  backscatter  mode,  it  is  interesting  to  note  that  at  98  and  140  GHz  the 
bistatic  scatter  signal  shows  an  angular  dependence  distinctly  more  influenced  by 
the  specular  mechanism  than  the  backscatter  signal.  This  can  be  qualitatively 
understood  in  the  following  way.  The  individual  scatterers  in  the  top  layer  of  the 
snow,  which  is  most  active  in  scattering,  have  certain  random  excursions  in  height 
from  the  average  surface  plane.  In  terms  of  the  Rayleigh  criterion  the  sum  signal 
from  the  individual  sources  will  be  more  coherent  the  lower  the  frequency  and  the 
lower  the  scattering  angle.  This  means,  for  a  given  snow  surface  and  frequency, 
bistatic  returns  are  more  coherent  as  seen  in  Figures  14  and  15.  In  the  extreme 

case  of  very  low  scattering  angles  this  has  been  demonstrated  bv  our  multipath 

12 

propagation  measurements. 

Transmission  measurements  through  this  snow  at  35  GHz  produced  a  result 
intermediate  to  that  measured  through  the  slab  of  sleet  and  the  slab  of  melt-freeze 
snow.  The  slab  containing  sleet  caused  almost  no  modulation  of  the  transmitted 
signal,  and  the  melt-freeze  snow  resulted  in  a  transmitted  signal  modulated  over 
the  entire  angular  range.  In  the  case  of  FT  metamorphic  snow  there  is  a  weak  in¬ 
crease  in  the  transmitted  signal  on  both  sides  of  perpendicular  incidence.  From  a 
grazing  angle  of  45  deg  on  down,  the  transmitted  signal  level  decreases  while  super¬ 
imposed  by  a  gentle  modulation.  Below  20  deg,  the  modulation  changes  to  a  series 
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of  abrupt  and  deep  fades.  The  increase  in  signal  level  at  verv  low  grazing  angle 
d.  e  to  diffraction  around  the  edge  of  the  basket  does  not  show  in  this  graph 
as  it  did  in  the  other  attenuation  measurements  at  this  and  the  higher  fre¬ 
quencies.  The  98-  and  140-GHz  plots  also  show  notable  differences  with  regard 
to  the  earlier  figures.  At  98  GHz,  attenuation  is  generally  lower  than  in  the  pre¬ 
vious  cases  with  the  exception  of  a  few  deep  nulls,  one  peculiarly  located  at  around 
TO  deg  grazing  angle.  The  fluctuation  rates  compare  somewhat  with  ‘hose  seen  for 
sleet  and  mclt-frceze  snow.  This  is  in  contrast  with  the  140-GHz  d.-  .  They  shew 

higher  attenuation  than  m  Figure  12  with  a  similar  number  of  fluctuations,  but 
similar  attenuation  as  in  Figure  la  with  a  lower  number  el  fluctuations.  It  is  appar¬ 
ent  from  this  discussion  of  sample  cases  that  results  are  not  always  systematic  or 
at  least  not  easily  explainable  in  terms  of  snow  parameters  observed. 

The  data  collected  after  the  bli/./ard  constitute  the  largest  number  of  cases 
with  only  limited  changes  in  snow  characteristics.  They  will  be  used  here  for  com¬ 
parisons  of  fluctuation  rates  and  signal  level  change  with  grazing  angle  as  pievious- 
lv  developed  from  the  models.  Based  on  the  appropriate  curve  of  Figure  10,  value 
for  the  number  of  backscatter  signal  fluctuations  per  deg  was  calculated  at  35,  98, 
and  140  GHz  at  grazing  angles  of  90,  45,  and  15  deg.  The  results  are  plotted  in 
Figure  If  as  three  dashed  curves,  each  representing  one  of  the  three  frequencies. 

In  the  same  figure,  the  three  solid  curves  represent  the  corresponding  fluctuation 
rates  derived  from  the  measured  data.  The  measured  values  were  arrived  at  bv 
counting  the  number  of  fluctuations  occurring  within  5  deg  of  either  side  of  grazing 
angles  of  90,  45,  and  15  deg.  Average  rates  were  those  derived  from  20  data  runs 
at  '35  GHz,  18  runs  at  98  Gllz,  and  14  runs  at  140  GHz,  One  observation,  when  in¬ 
specting  Figure  17,  is  a  measured  fluctuation  rate  greater  than  that  predicted  at 
all  frequencies  and  at  all  grazing  angles.  This  is  particularly  evident  at  the  lower 
grazing  angles  where,  for  example,  the  ratio  of  measured  and  predated  rates 
approaches  4  at  140  GHz  and  15  deg.  The  model  used  to  derive  backscatter  fluc¬ 
tuation  rates  assumes  single  scattering,  it  also  considers  unlv  the  phase  interfer¬ 
ence  between  a  scatterer  at  a  fixed  distance  d  from  the  center  of  rotation  and  a 
seatterer  immobile  at  the  center  of  rotation.  This  would  tend  to  represent  the  true 
ease  at  perpendicular  orientation  of  snow  surface  and  transmitter  beam.  When  the 
slab  tilts  relative  to  the  beam  for  a  lower  grazing  angle,  the  edge  of  the  3-dH  beam 
corresponds  to  a  larger  d  (the  larger  half-axis  of  the  elliptical  antenna  footprint). 

At  a  grazing  angle  of  15  deg,  d  has  increased  by  a  factor  of  2.  7  in  this  direction. 

If  we  further  consider  that  the  radar  range-factor  puts  the  closest  portion  of  the 
footprint  at  a  noticeable  advantage  in  terms  of  power  received  per  unit  scattering 
area,  this  could  account  for  the  discrepancy  between  calculated  and  measured  fluc¬ 
tuation  rates.  A  more  elaborate  model  is  in  order  for  the  backscatter  case.  Mul¬ 
tiple  scattering  might  be  a  contributing  factor,  but  is  not  required  to  explain  the 
experimental  results.  Averaging  of  the  number  of  fluctuations  within. i  10-deg 
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interval  is  a  subjective  procedure, 
observed. 


It  must  account  for  part  of  the  differences 
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Because  of  tile  more  limited  angular  range,  over  which  bistatic  scatter  data 
were  observable,  no  attempt  was  made  to  verify  the  model  fluctuation  rates  in 
figure  10.  It  is  clear  in  the  bistatic  case  that  the  elliptical  shape  of  the  scattering 
area  on  snow  also  leads  to  higher  fluctuation  rates  at  low  grazing  angles  than  cal¬ 
culated  for  a  circular  area. 

It  has  been  demonstrated  bv  the  sample  cases  that  upon  transmission  through 
snow,  the  received  signal  level  varies  in  a  complex  manner  with  changing  direction 
of  penetration.  We  had  modeled  four  different  mechanisms  in  Figure  10  which 
could  explain  interaction  of  multiple  signals  in  the  transmission  mode.  At  low- 
grazing  angles,  a  signal  diffracted  around  the  outer  edge  of  the  basket  compares 
well  with  that  transmitted  at  an  oblique  angle  through  the  snow.  The  oblique  path 
through  the  snow  is  longer  than  the  thickness  of  the  slab.  It  depends  on  the  type  of 
snow,  thickness,  and  frequency  at  which  basket  angle  the  diffracted  signal  becomes 
a  significant  portion  of  the  total  signal  received.  Diffraction  effects  must  be  seen 
as  an  undesired  side  effect  of  our  measurement  setup.  For  each  plot  of  the  post¬ 
blizzard  data,  a  grazing  angle  was  established  below  which  the  fluctuations  were 
attributed  to  diffraction.  The  number  of  fluctuations  was  counted.  Average  numbers 


for  all  attenuation  plots  from  this  period  are  listed  in  Table  4.  Measured  numbers 
of  fluctuations  are  those  between  10  deg  and  the  average  angle  of  onset.  Note  that 
this  latter  angle  increases  with  frequency  as  can  be  expected.  Predicted  numbers 
are  based  on  an  average  rate  of  range  change  of  1.  9  mm/deg,  estimated  from  fig¬ 
ure  10  in  the  diffraction  region.  Remember  that  the  curve  labeled  Reflection'' 
approximates  the  one  for  diffraction.  Specific  numbers  in  Table  4  represent 
1.9  times  the  deg-range  of  diffraction,  divided  by  wavelength.  The  results  sup¬ 
port  the  assumption  of  diffracted  signals. 


Table  4.  Fluctuation  Characteristics  of  Attenuated  Signal 


Frequency 

(GHz) 

35 

98 

140 

Diffraction 

Region 

Average  Angle  of 

Onset  (deg) 

22.  2 

27.  3 

29.  4 

Measured  Number 
of  Fluctuations 

3.  1 

13.  3 

18.  4 

Predicted  Number 
of  Fluctuations 

2.  7 

11.0 

18.  5 

30-  to  90-  to  60-deg 
Region 

Measured  Number  of 
Slow  Fluctuations 

15.  9 

24.  0 

Predicted  Number  of 
Slow  Fluctuations 

2.  7 

7.  8 

11.8 

30-  to  90-  to  60-deg 
Region 

Measured  Number  of 
Fast  Fluctuations 

54.  0 

84.  2 

Predicted  Number  of 

Fast  Fluctuations 

_ 

20.  9 

60.  0 

90.  0 

The  other  three  mechanisms,  suggested  for  the  explanation  of  observed  fluctua¬ 
tions,  are  single  and  multiple  scattering  and  reflection  off  the  snow/air  interfaces, 
if  one  mechanism  were  predominant,  measurements  should  identify  it  bv  its  charac¬ 
teristic  angular  dependence.  Unfortunately,  results  do  not  bear  this  out  very  well. 
If  signals  due  to  scattering  and  reflection  reach  comparable  amplitudes,  the  angular 
variation  becomes  inconclusive.  In  Table  4  we  have  therefore  shown  total  numbers 
of  fluctuations  over  the  90-deg  interval  of  basket  rotation  outside  the  diffraction 
region.  Measured  data  are  again  averages  over  all  available  plot  s.  As  it  gener¬ 
ally  appeared  that  fluctuations  could  be  grouped  into  fast  and  slow  ones,  these  have 
been  counted  separately.  It  is  not  clear  that  single  and  multiple  scattering  should 
lead  to  such  a  distinctive  distribution  of  fluctuation  rates.  For  the  sake  of  a  simple 
comparison,  we  have  assumed  that  "Single  Scatter"  in  the  90-  to  90-  to  60-deg 


region  predicts  an  upper  bound  of  0.  25  nim/deg  for  the  rate  of  range  change.  The 
corresponding  number  of  "Multiple  Scatter"  is  2.  1  mm/deg.  Multiplying  these 
numbers  by  the  deg-range  of  observation  and  dividing  bv  wavelength  yields  the  pre¬ 
dicted  numbers  of  slow  and  fast  fluctuations,  respectively. 

Table  4  omits  measured  data  in  the  30-  to  90-  to  60-deg  region  at  35  Cillz. 
Fluctuations  were  non-existing  or  too  inconsistent  to  arrive  at  meaningful  values 
at  this  frequency.  At  98  GHz  and  140  GHz  the  number  of  observed  fast  fluctuations 
is  indeed  bounded  by  the  prediction.  However,  observed  numbers  of  slow  fluctua¬ 
tions  are  higher  by  a  factor  of  2  or  more  than  predicted.  No  further  conclusions 
will  be  drawn  from  this,  since  the  model  apparently  does  not  match  the  complexity 
of  the  real  situation. 

4.2  Summary  of  Experimental  Data 

While  discussing  the  results  of  all  scatter  and  attenuation  measurements  ob¬ 
tained  on  snow  samples,  reference  will  be  made  to  meteorological  observations 
and  measurements  of  physical  snow  parameters  or  ground  truth  data  as  they  are 
culled  in  reference  to  millimeter-wave  imaging  systems.  These  data  are  compiled 
into  tables.  In  the  appendix,  one  table  is  given  for  each  of  the  20  days,  on  which 
scatter  and  attenuation  measurements  were  made.  The  data  underline  the  struc¬ 
tural  and  developmental  variability  of  a  medium  which  is  not  obvious  to  the  casual 
observer.  The  surface  of  snow  cover  may  also  vary  substantially  in  terms  of  the 
wavelengths  of  interest  here.  For  illustration  purposes,  two  of  the  slabs  used  in 
the  experiments  are  shown  in  Figures  18  and  19.  The  surface  of  a  slab  cut  on 
2  March  from  FT  snow  is  almost  smooth  (Figure  18).  In  contrast.  Figure  19  dis¬ 
plays  a  slab  with  a  heavily  cusped  surface.  It  was  cut  on  2  February  and  owes  its 
surface  structure  to  MF  metamorphism. 

All  data  presented  in  this  section  are  average  peak  values.  Referring  back  to 
the  original  data  plots  in  Figures  12  through  15,  the  only  practical  wav  to  manually 
derive  average  signal  levels  from  these  piots  was  to  average  the  local  peaks  of  the 
rapidly  fluctuating  signals.  An  analysis  of  how  these  average  peak  signals  relate 

to  mean  signal  levels  will  be  included  in  the  forthcoming  report  on  in-situ  snow 
14 

scatter  measurements.  Hear  in  mind  that  all  subsequent  graphs  in  this  report 
pertaining  to  actually  measured  data  represent  averages  of  peaks. 

Based  on  the  mode  of  operation  of  the  equipment  (continuous  measurements  are 
made  while  the  basket  rotates),  data  exist  for  a  contiguous  range  of  grazing  angles 
Q .  In  the  in-situ  configuration,  grazing  angles  had  to  be  chosen  parametrically  to 
limit  the  total  number  of  measurements.  Values  ford  selected  were  90,  45,  and 
15  deg.  The  lowest  angle  over  terrain  is  dictated  both  by  system  sensitivity  and 
the  physical  limits  of  the  snow  field. 
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Figure  20.  Backscatter 
Coefficient  From  Sleet, 
(data  of  16  January  1978, 
vertical  polarization) 
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On  17  January,  the  condition  of  the  snow  was  essentially  unchanged  from  that 
of  the  previous  day.  The  backscatter  experiment  was  repeated,  again  at  vortical 
polarization.  Measurements  of  a°  at  45  and  15  deg  were  almost  unchanged.  At 
90  deg,  there  was  a  5-dB  spread.  At  35  GHz,  <j°  was  found  to  be  11  dll,  hardly 
changed  from  the  previous  day.  At  98  GHz,  a  decrease  bv  3  dB  to  8  dB  was  ob¬ 
served.  At  140  GHz  the  change  was  positive  by  2.  5  dB  to  13.  5  dB. 

On  18  January,  a  heavy  rain/snow  storm  took  place  which  deposited  5  to  8  cm 
of  wet,  irregularly  shaped  ice  grains  on  top  of  the  existing  snow  cover.  The  mean 
diameter  of  the  grains  was  approximately  1  mm.  Following  this,  a  severe  snow¬ 
storm  hit  on  20  January.  A  record  snow  accumulation  for  a  24-hr  period  of  56  cm 
was  measured  in  Boston.  No  rain  was  associated  with  this  storm  which  was  followed 
by  clear,  cold  (subfreezing)  days.  The  next  set  of  oackseatter  measurements  were 
performed  on  24  January.  The  backscatter  coefficients  measured  on  this  day  are 
shown  in  Figure  21.  The  snow  parameters  of  the  slab  used  are  given  in  the  appen¬ 
dix.  The  sample  had  been  cut  for  use  on  19  January.  However,  it  was  not  used  as 
tests  were  stopped  in  order  to  modify  the  basket  drive  mechanism.  I'ntil  testing 
resumed  on  21  January,  the  slab  lay  horizontally  on  a  platform  at  the  test  site. 

During  this  time,  snow  from  the  20  January  storm  was  deposited  on  the  slab.  A 
top  layer  formed,  3.  5-em  thick,  with  a  density  of  approximately  0.25  g/cm  *.  t'nder 
the  new  snow,  a  3-cm  layer  of  irregularly  shaped  ice  grains  of  1  to  2  mm  diameter 
had  a  density  of  0.  38  g/em'  .  This  layer  had  been  deposited  in  the  rain/snow  storm 
of  18  January.  Fnder  this  lay  a  4.  5-em  thick  layer  of  the  dense-  (0.53  g/etn  S  snow, 
used  in  the  tests  of  16  and  17  January. 
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Figure  21.  Backscatter 
Coefficient  From  Melt-Freeze 
Snow,  (data  of  24  January  1978, 
vertical  polarization) 


There  is  a  marked  change  in  the  angular  dependence  of  the  backscatter  coeffi¬ 
cient  especially  at  35  GHz.  At  this  frequency,  at  a  grazing  angle  of  90  deg,  a0  has 
dropped  bv  5  dB  relative  to  the  level  of  16  January.  At  the  shallower  angles  of  45 
and  15  deg,  the  backscatter  coefficient  has  increased  by  6  and  3  dB,  respectively. 
The  newer,  less  dense  top  snow  had  removed  all  indications  of  a  specular  nature 
of  the  backscatter  coefficient.  There  is  a  lesser  change  in  o°  at  98  and  140  GHz. 

At  90  deg,  the  value  of  (7°  is  essentially  unchanged.  At  45  and  15  deg,  the  98-GHz 
value  has  increased  bv  approximately  3  dB.  At  140  GHz,  a°  has  ircreased  by 
3  dB  at  15  deg  but  is  unchanged  at  45  deg. 

Main  fell  on  25  and  26  January,  followed  by  sunny  and  cold  days  up  to  6  Feb¬ 
ruary  when  the  tests  were  stopped  due  to  heavy  snowfall.  Additional  backscatter 
measurements  were  made  on  27,  28,  30,  31  January  and  1,  2,  4,  6  February.  The 
sample  slabs  were  cut  from  the  upper  10  cm  of  snow  cover  and,  with  the  exception 
of  a  lower  ice  layer,  consisted  entirely  of  snow  deposited  in  the  storm  of  20  January. 
To  determine  the  effect  of  the  ice  layer,  some  of  the  tests  were  repeated  following 
removal  of  the  ice  layer.  No  change  in  O'  was  detected. 

According  to  the  tables  in  the  appendix,  the  volume  characteristics  (density, 
temperature,  and  so  on)  of  the  snow  did  not  vary  larch  during  this  period.  The 
rainstorms  of  25  and  26  January  did  cause  cusp-shu,.<  d  depressions  to  be  formed 
on  the  snow  surface.  An  example  of  this  phenomeno"  was  shown  in  Figure  19.  The 
cusps  were  about  1  cm  deep  and  about  5  cm  in  diameter.  Some  11  cusps  were 
counted  per  meter.  Nonetheless,  the  snow  return  did  not  change  appreciably  during 
this  period.  This  is  illustrated  in  Figure  22.  The  average  of  the  data  collected  on 
16  and  17  January  is  shown  as  solid  lines.  The  dashed  lines  are  the  result  of 
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averaging  all  data  collected  in  the  period  from  24  January  to  t>  February.  Com¬ 
paring  the  dashed  curves  in  Figure  22  with  the  data  of  24  January  in  Figure  21, 
one  sees  only  small  changes  at  each  grazing  angle  and  frequency.  In  contrast,  the 
differences  between  the  solid  and  dashed  curves  in  Figure  22  are  more  significant. 
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Figure  22.  Average  Hack- 
scatter  Coefficient  From  Sleet 
( —  data  of  16  and  17  January 
1978)  and  From  Melt-1  ree/e 

Snow,  ( - data  of  24  January 

to  6  February  1978),  Vertical 
Polarization 
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The  test  program  was  expanded  following  the  blizzard  of  t>  February.  In  addi¬ 
tion  to  collecting  data  at  vertical  polarization,  measurements  were  made  at  hori¬ 
zontal  polarization  and  in  the  cross-polarized  mode.  During  the  two- week  period 
from  the  end  of  the  storm  to  the  resumption  of  the  measurements  on  21  February, 
metamorphism  proceeded  primarily  by  the  FT  route.  This  is  a  slow  process  com¬ 
pared  to  that  brought  about  by  the  heavy  rain/ snow  storms  in  January.  However, 
it  appears  that  two  weeks  were  sufficient  for  the  process  to  reach  an  equilibrium 
state. 

With  the  exception  of  a  test  to  be  discussed  later,  all  tests  conducted  in  this 
period  used  snow  samples  which  wore  12.7  cm  thick.  Keference  to  the  ground  truth 
data  (see  Appendix  A)  shows  that  the  top  2.5  cm  of  each  sample  consisted  of  a  layer 
of  crust.  This  layer  was  produced  by  solar  heating  and  melting  of  the  snow  during 
the  day  and  subsequent  refreezing  at  night.  The  remaining  snow  below  this  layer 
was  composed  of  grains  produced  by  FT  metamorphism.  The  grains  were  only 
weakly  bound  together,  producing  a  fragile  sample  which  was  supported  mostly  bv 
the  crustal  layer.  The  density  in  these  samples  (0.  26  g/cm'S  was  roughly  the  same 
as  in  the  samples  used  in  January.  The  grain  sizes  were  also  similar. 
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Data  wore  gathered  on  7  days  in  the  period  from  21  February  through  2  March. 
Faeh  day  the  sequence  of  scatter  and  transmission  measurements  at  polarizations 
V  Y .  Hll,  and  HY  started  with  a  different  frequency.  The  sequence  was  then  re¬ 
peated  at  the  other  two  frequencies.  A  single  snow  sample  was  used  for  each  day  's 
tests.  Following  the  completion  of  the  measurements,  .<  run  duplicating  the  initial 
one  of  the  day  was  made.  The  purpose  was  to  establish  that  no  ma  jor  change  in  the 
sample  properties  had  occurred  during  the  test  period.  A  different  frequency  was 
chosen  on  a  rotational  basis  to  begin  the  tests  on  different  days. 

Since  tile  properties  of  the  samples  did  not  vary  much  from  day  to  day,  data 
gathered  on  individual  days  will  not  be  considered  separately.  Instead,  like  data 
from  each  day  have  been  averaged  together.  The  average  backscattor  coefficient 
as  a  function  of  the  grazing  angle  is  presented  in  Figure  23  .or  each  frequency  and 
polarization.  Compare  the  backscatter  coefficient  for  vertical  polarization  (dashed 
lines  in  Figure  23)  with  the  corresponding  data  in  the  January  tests  (dashed  lines 
in  Figure  22).  The  35-til Iz  data  in  each  figure  show  the  backscatter  coefficients  at 
JO  deg  as  equal.  At  -15  and  15  deg,  average  backscatter  returns  measured  from 
the  February  snow  samples  at  35  till/  are  3  to  1  dll  below  those  measured  from  the 
January  samples.  At  08  till/.,  t lie  average  return  from  the  February  snow  is 
approximately  3  dll  lower  al  each  grazing  angle.  The  return  at  140  till/,  off  Feb¬ 
ruary  samples  is  1  to  2  dll  higher  than  that  measured  in  January. 


a  45 
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Figure  23.  Average  Back¬ 
scatter  Coefficient  From 
Fquitemperature  Snow,  (data 
of  21  February  to  2  March 
1978,  polarizations  Hll  — , 

V A  ---,  IlY  -  •  -  ■  ) 


The  average  35-CiIlz  backscatter  coefficient  at  horizontal  polarization  exceeds 
the  corresponding  values  at  vertical  polarization  by  about  2  dll.  At  98  GHz,  the 
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values  of  these  quantities  are  roughly  equal.  The  average  values  of  the  backscatter 
coefficient  for  horizontal  polarization  at  140  GHz  is  approximately  2  dB  less  than 
that  for  vertical  polarization  at  all  grazing  angles.  The  cross-polarized  back¬ 
scatter  coefficient  (J4^,  displays  an  interesting  behavior.  At  35  GHz,  a°^y  is 
18  dB  lower  than  at  a  grazing  angle  of  90  deg  and  12  dB  lower  at  45  and 

15  deg.  The  angular  dependence  of  and  °°y\j  at  35  GHz  is  somewhat  specular 

(similar  to  the  early  sleet  samples),  and  the  small  levels  of  the  cross-polarized 
component  relative  to  the  like-polarized  levels  are  further  indication  of  scattering 
from  a  surface  which  appears  smooth  at  the  particular  wavelength  (35  GHz).  A 
dramatic  change  in  is  seen  as  the  frequency  is  increased.  The  strength  of 

the  cross-polarized  return  relative  to  the  like -polarized  return  goes  up  with  fre¬ 
quency.  At  98  GHz,  the  cross-polarized  a°  is,  on  the  average,  only  5  dB  below 
the  like-polarized  <x°.  At  the  highest  frequency,  140  GHz,  the  average  cross- 
polarized  and  horizontal  like-polarized  backscatter  coefficients  are  almost  equal. 

As  has  already  been  mentioned,  the  snow  samples  used  in  the  February  tests 
were  covered  by  a  2-  to  3-cm  thick  crust.  To  determine  the  strength  of  the  back¬ 
scatter  from  the  crust  relative  to  that  from  the  full  slab,  measurements  were  made 
in  one  case  on  the  crust  alone.  The  sample  was  cut  to  conform  to  the  standard 
thickness  for  this  series  (12.  7  cm)  and  the  backscatter  coefficient  determined  at 
98  GHz.  Following  this  run,  the  bottom  10-cm  of  fragile  snow  was  removed  from 
the  sample,  leaving  only  the  crust  in  place.  Then,  the  backscatter  measurement 
was  repeated  using  the  thinner  sample.  After  this,  the  thin  sample  was  used  to 
measure  the  backscatter  coefficients  at  35  and  140  Gllz.  The  results  at  98  GHz 
using  the  original  12.  7-cm  sample  agreed  to  within  1  dB  with  those  measured  on 
the  thin  (2.5  cm)  sample.  Since  98  (ill/  was  the  only  frequency  at  which  a  direct 
comparison  was  possible,  at  the  other  two  frequencies  a  comparison  was  made 
against  averages  of  samples  consisting  of  ,  rust  and  underlying  fragile  snow  from 
the  whole  period.  To  illustrate  the  general  ♦rend  i.ln..  of  o ^  ,,nd  ob¬ 

tained  from  the  crust  only  have  been  .\ .  r  .ged  t.  g.'n.  i.  in  I  igurt  24,  these  val¬ 
ues  connected  by  dashed  lines  art  plotted  to:  „  ,-n  .  is  .  function  of  the 

grazing  angle.  The  averages  of  o and  ’  I  ,k  .  ■  2  ■  Iron,  the  thicker 

samples  are  shown  in  Figure  24  b\  solid  lira  dl  ’tuna  rust  comparison 

against  the  average  of  the  thiek-samp).  a  >t  .  .  ,  •  •  mi  nt  like  that  against 

the  single  thick  sample.  This  suggests  tha-  ,  .  *♦,  ting  is  not  taking 

place  at  depths  beyond  2.  5  cm.  Attenuation  .  >.t  !!t  >.  n*  ’in  n«  vt  sta  tion  support 
the  assumption.  At  35  GHz,  deeper  penetration  into  th-  -1  .8  mast  bo  exported,  and 
the  higher  backscatter  coefficients  shown  in  t  igurt  24  air  plausible.  There  is  no 
simple  explanation  as  to  why  the  140-GH/  results  should  show  the  same  trend.  This 
is  particularly  so,  since  the  98-Gllz  results  indicate  that  the  effect  of  the  fragile 


snow  was  small.  It  must  be  kept  in  mind  that  this  was  only  a  single  case  of  thin 
snow  crust. 


Figure  24.  Comparison  of 
Baekseatter  Coefficient 
Measured  off  Thick  and 
Thin  Snow  Slabs, 

(12.7-cm  slab  — , 

2.  5-cm  slab  -  -  - ) 
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The  angular  dependence  of  scattered  signal  intensity  as  a  function  of  grazing 
angle  according  to  Clapp  s  models1^  could  shed  additional  light  on  the  scattering 
mechanisms  of  metamorphic  snow.  A  volume-scattering  medium  would  be  dis¬ 
tinguishable  by  an  angular  a  "-dependence  following  a  sin  6  -law  .  A  surface- 

*> 

scattering  medium  would  display  a  sin“  $  -dependence  on  grazing  angle.  Specular 
reflection  would  be  indicated  by  an  angular  dependence  proportional  to  the  product 
of  the  antenna  patterns.  Referring  back  to  Figure  11,  a  volume-scattering  medium 
fitting  Clapp's  model  exhibits  a  decrease  in  scattered  power  of  1.  B  dB  between 
grazing  angles  of  90  and  45  deg.  At  15  deg,  the  scattered  power  has  dropped  to 
7. 3  dB  below  that  for  perpendicular  incidence.  Hough-surface  scattering  according 
to  Clapp  s  model  is  indicated  by  a  signal  drop-off  of  2.  9  and  11.5  dB  at  the  lower 
grazing  angles.  True  specular  scattering  measured  with  the  average  beam  width 
used  in  the  scatterometers  requires  a  signal  reduction  bv  10  dB  within  t>  deg  from 
perpendicular  incidence. 

Theoretical  curves  were  plotted  and  overlaid  onto  the  average  baekseatter 
curves  for  sleet  and  MF  snow  in  Figure  22  and  for  FT  snow  in  Figure  23.  The 
three-point  representation  of  angular  dependence  of  the  average  baekseatter  coeffi¬ 
cient  did  not  preserve  the  real  drop-off  observed  on  sleet  as,  by  way  of  example, 
in  Figure  12.  There  it  was  shown  that  a  true  specular  response  is  indeed  possible 
at  35  GHz.  The  value  of  <7°  at  45  deg  is  at  least  IB  dB  below  that  for  nonspecula r 


scatter,  ruling  out  the  latter  mechanism.  We  saw  in  the  higher-frequency  signals 
from  sleet  in  Figure  12  a  much  stronger  incoherent  component  in  the  received 
signal.  Comparing  the  average  levels  in  Figure  22  with  the  models,  there  exists 
still  a  distinctive  difference  of  8  dll  or  more  at  a  45-deg  grazing  angle,  attesting 
to  a  strong  specular  component  in  the  scattered  signal.  The  melt-freeze  snow 
averages  in  the  same  figure  (the  dashed  lines)  conform  in  their  angular  dependence 
much  more  closely  to  the  incoherent  scattering  models.  One  may  have  to  exempt 
the  35-Gllz  curve  from  this  statement,  since  between  90  deg  and  45  deg  the  slope 
of  the  curve  is  still  steeper  than  expected  from  the  model.  Even  though  the  sample 
in  Figure  13  does  not  show  a  distinctive  "main  lobe",  there  is  on  the  average  a 
considerable  contribution  from  a  specular  component  at  this  frequency  and  over 
this  type  of  snow.  Between  45  and  15  deg,  the  slope  of  the  35-GHz  line  agrees  well 
with  the  one  calculated  for  volume  scatter.  Fitting  the  model  curves  to  the  mea¬ 
sured  data  at  the  upper  two  frequencies  yields  good  agreement  in  the  angular  range 
between  90  and  45  deg.  That  is,  any  indication  of  specular  scatter  is  clearly  ruled 
out  in  this  angular  range.  The  differences  between  expected  values  for  volume 
scatter  and  rough-surface  scatter  are  not  significant  enough  to  conclude  which 
scatter  mechanism  is  predominant  above  a  grazing  angle  of  45  deg.  Below  45  deg, 
and  judging  by  the  slope  of  the  lines  connecting  the  15-  and  45-deg  points,  volume 
scatter  is  the  mechanism  accounting  for  a°.  Of  course,  the  lowest  grazing  angle 
considered  here  is  not  low  enough  to  produce  a  really  convincing  result.  It  is  well 
possible  that  contributions  from  both  incoherent  scatter  mechanisms  exist  with 
comparable  magnitude,  making  it  again  difficult  to  sort  out  one  mechanism  or  the 
other. 

Figure  23,  showing  the  results  for  ET  snow,  is  not  too  different  in  angular 
dependence  of  backscatter  coefficient  from  Figure  22.  The  35-GHz  data  show  a 
behavior  intermediate  to  sleet  and  MF  snow,  presumably  because  the  crustal  layer 
approaches  the  mirror  properties  of  the  frozen  sleet.  At  angles  below  45  deg,  the 
line  slopes  down  in  agreement  with  the  one  postulated  for  volume  scatter.  The  same 
is  true  at  98  and  140  GHz.  The  crossover  found  at  98  GHz  between  VV  and  1IH 
curves  most  likely  is  not  significant.  The  two  different  slopes  that  occur  are  still 
best  identified  with  that  of  volume  scatter.  Between  45  and  90  deg,  slopes  measured 
agree  with  the  calculation. 

4.2.2  T  HANSM1SSION  DATA 

In  Section  4.  1  the  nature  of  the  received  signal  was  discussed  when  the  system 
operated  in  the  transmission  mode.  Transmission  data  will  be  used  to  obtain  a 
measure  of  snow  attenuation  at  millimeter  wavelengths.  The  analysis  will  be 
limited  to  determining  an  attenuation  coefficient  for  perpendicular  incidence.  Tor 
35  GHz,  this  quantit”  is  derived  directly  from  the  transmitted  signal  at  perpendicular 


incidence.  In  the  case  of  the  98-  and  140-GHz  data,  where  the  signal  fluctuates 
rapidly,  the  average  peak  transmitted  signal  was  determined  by  averaging  the 
transmitted- signal  peaks  measured  within  5  deg  on  either  side  of  perpendicular 
incidence.  This  corresponds  to  the  procedure  used  to  determine  the  average  peak 
backscatter  as  a  function  of  grazing  angle. 

All  attenuation  data  from  the  period  before  the  6  February  storm  are  grouped 
together  in  Table  5,  and  all  data  collected  subsequent  to  that  day  in  Table  6.  Ir. 
each  table  the  date,  the  snow  thickness,  and  the  density  are  given.  One  column 
shows  the  average  attenuation  measured  at  each  frequency.  In  the  last  column  the 
attenuation  coefficient  expressed  ir.  dll  per  meter  is  given.  All  measurements 
shown  in  Table  5  were  made  in  the  vertical  polarization  mode,  following  the 
blizzard,  the  attenuation  was  measured  in  both  polarization  modes  and  both  are 
shown  in  fable  8.  Under  field  conditions  it  is  reasonable  to  assume  a  measurement 
accuracy  for  attenuation  no  better  than  2  dii.  Attenuations  below  this  level  were 
not  used  to  obtain  average  values.  A  fair  amount  of  the  35-Gllz  data  fall  into  this 
category.  They  are  shown  in  parentheses. 

To  provide  a  basis  for  comparison,  attenuation  coefficients  have  been  calcu¬ 
lated  assuming  snow  to  be  a  dielectric  medium  composed  of  ice  and  air.  Miles 
1 B 

and  Uluby  give  the  attenuation  coefficient  through  snow  as 

8.  68  ,7  i  . 

I.  - - dU/m  ,  1. 14) 

with  t  as  the  real  part  and  t  .  the  iniaginarv  part  of  the  dielectric  constant  c  of 
t  i  ( '■>  -•> 

the  mixture.  In  the  snow  multipath  report  "  we  derived  the  dielectric  constant  of 

snow  from  a  mixing  formula  cited  by  Gumming.  Stiles  and  L'i-bv  use  Wiener  s“° 

21 

mixing  formula  referenced  by  Evans  which  permits  one  to  take  internal  structure 
of  snow  into  consideration.  If  we  assume  the  dielectric  constant  of  air  to  be  unity, 
Wiener's  equation  reduces  to 

(i  -  1)  <t  f  f)  =  f  (t  -  1 )  (t  +  F)  ,  (15) 

s  s 
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<  dielectric  constant  of  ice, 

f  fractional  volume  occupied  by  ice, 

!■'  form  factor,  ranging  from  zero  to  infinity; 

F  =  2  for  ice  spheres. 

The  dielectric  constant  of  ice,  a  low-loss  medium,  approximately  relates  to  the 
real  and  imaginary  component  (n^,  n.)  of  the  refractive  index  given  bv  Kay  1  and 
quoted  earlier  in  the  report  as 

■> 

t  -  n  ~  +  i  2  n  n .  .  (10) 

r  '  r  l 

Eased  on  a  form  factor  F  2,  dielectric  constants  of  snow  have  been  calculated  from 
\\  iener  s  formula  at  the  three  frequencies  for  average  densities  measured  during 
the  tests.  Attenuation  coefficients  for  the  various  snows  and  for  plain  ice  were  then 
derived  using  Eq.  (14).  The  results  are  tabulated  m  Table  7.  A  comparison  of 
experimental  data  in  Tables  5  and  6  with  the  calculated  ones  in  Table  7  makes  clear 
that  millimeter- wave  attenuation  through  snow  which  apparently  involves  particle 
resonances  is  not  described  satisfactorily  by  the  Wiener  formula.  Even  when  im¬ 
plementing  the  equation  with  F  equal  to  infinity,  the  predicted  attenuation  coeffi¬ 
cients  come  nowhere  near  the  ones  measured  at  98  and  140  till/..  The  same  can  be 
said  w'hen  applying  the  Gumming  formula. 

Examining  Tables  5  and  ti  in  more  detail,  one  sees  quite  a  bit  of  scatter  m  the 
data.  However,  a  number  of  trends  are  discernible.  First,  the  attenuation  coeffi¬ 
cient  measured  at  35  Gil/,  through  the  post-blizzard  snow  (Table  6,  ET  snow  )  is 
intermediate  to  that  measured  through  sleet  and  1UF  snow  (Table  5).  Second,  this 
does  not  appear  to  be  the  case  at  the  two  higher  frequencies.  With  the  exception  of 
24  February,  the  attenuation  coefficients  measured  at  08  and  140  GHz  through  ITT 
snow  are,  on  the  whole,  lower  than  all  values  obtained  from  sleet  and  WF  snow. 
Third,  it  appears  that  the  attenuation  coefficient  may  depend  on  the  slab  thickness. 

In  the  pre-blizzard  tests  of  Table  5,  the  largest  values  of  1.  were  measured  at  08 
and  140  GHz  through  the  sleet  slab  of  17  January.  It  was  also  the  thinnest  of  all  the 
slabs  tested  during  this  period.  At  35  GHz,  the  attenuation  t  hrough  this  slab  wuis 
too  small  to  be  considered.  Also,  the  98-  and  140-GHz  attenuation  coefficients 
through  the  6-cm  slabs  seem  higher  on  the  average  than  those  through  10  to  14-oni 
slabs.  This  conclusion  is  strengthened  by  attenuation  measurements  performed  on 
24  February  on  the  single  2.  5-cm  crust  slab  vs  the  others  during  this  period  of 
stable  snow  structure  which  had  12.  7  cm  thickness.  The  particle  size  and  density 
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contained  in  the  crust  did  not  differ  appreciably  from  those  in  the  underlayer.  It 
is  not  clear,  if  the  differences  in  hardness  between  the  top  and  lower  layers  suffice 
to  make  them  different  in  terms  of  quantity  !..  The  brittle  portion  was  not  investi¬ 
gated  separately  for  its  attenuation  coefficient. 

The  theory  of  multiple  scattering  in  optics'"  provides  a  reasoning  for  In¬ 
dependence  on  sample  thickness.  Consider  a  beam  of  photons  incident  upon  a 
scattering  medium.  As  the  beam  propagates  into  the  medium,  photons  are  scat¬ 
tered  out  of  it  resulting  in  attenuation  of  the  beam.  If  the  density  of  scatterers  is 
high  enough  and  the  medium  is  thick  enough,  eventually  some  of  the  photons  scattered 
out  of  the  beam  will  be  directed  back  in  through  multiple  scattering.  W  hen  this  is 
the  case,  the  net  rate  at  which  piiotons  arc  lost  out  of  the  beam  is  reduced.  Thus 
the  average  attenuation  coefficient  measured  from  a  sample  thick  enough  that  mul¬ 
tiple  scattering  is  appreciable  would  be  smaller  than  that  measured  through  a  thin 
slab  where  multiple  scattering  is  negligible.  A  sketch  of  this  behavior  appears  in 
Figure  25.  The  slope  of  the  curve  reduces  from  its  original  value  In  proportion 
to  the  effect  that  multiple  scattering  has  as  the  sample  thickness  increases.  Some 
steady  state  may  be  reached  when  the  amount  of  power  scattered  back  into  the  beam 
is  equal  to  that  scattered  out.  When  this  condition  is  reached,  the  slope  of  the 

Conceivably,  Wiener's  formula  is  a  closer 


Figure  25.  Hypothetical 
Attenuation  Curve  for 
Scattering  Medium 


curve  could  be  due  to  absorption  only, 
match  at  this  point. 


22.  Ishimaru,  A.  (1978)  Wave  Propagation  and  Scattering  in  Kandom  Media, 
Vol.  2,  Academic  Press,  New  York, 


Table  7.  Calculated  Attenuation  Coefficient  of  Snow  Using  Wiener  Formula 


Snow  Condition 

Density 

F  requency 

Dielectric  Constant 

Attenuation 

Coeff. 

(g/cm3) 

(GHz) 

e 

r 

£i 

(dB/m) 

Equitemperature 

0.  26 

35 

1.  37 

3. 5  X  10'* 

0.  96 

Snow 

98 

1.37 

2.  0  X  10-4 

1.56 

21  Feb.  -  2  March 

140 

1.  37 

1.  2  X  10‘4 

1.40 

Melt-F  reeze 

0.  29 

35 

1.42 

4.  1  X  10‘4 

1.  08 

Snow 

98 

1.42 

2.  3  X  10'4 

1.75 

21  Jan.  -  6  Feb. 

140 

1.42 

1.4  X  10'4 

1.  55 

Layered  Sleet 

0.40 

35 

1.  60 

6.  2  X  10~4 

1.  56 

and  Snow 

98 

1.60 

3.  5  X  10‘4 

2.  52 

24  January 

140 

1.  SO 

2.  1  X  10'4 

2.  28 

Sleet  Layer 

0.  53 

35 

1.86 

9.4  X  10'4 

2.  19 

17  January 

98 

1.86 

5. 3  X  10  ; 

3.53 

140 

1.  86 

3.  2  X  10  4 

3.  16 

Ice 

1.  00 

35 

3.  17 

3.  2  X  10‘3 

5.70 

98 

3.  17 

1.8  X  10~3 

9.  10 

140 

3.  17 

1.  1  X  10‘3 

8.  20 

In  line  with  the  foregoing  discussion,  a  least-mean-square  curve-fit  analysis 
of  the  attenuation  data  was  performed.  At  each  fi  quency  the  pre-blizzard  (melt- 
freeze)  and  post-blizzard  (equitemperature)  snow  were  considered  separately. 
Linear  and  quadratic  curve  fits  were  made  based  on  a  subjective  decision  as  to 
which  best  matched  the  data.  In  the  case  of  35  GHz,  straight-line  fits  proved  to  be 
best.  A  parabolic  fit  was  selected  at  the  higher  frequencies.  For  all  data  sets, 
zero  attenuation  for  zero  slab  thickness  was  included  as  one  input  data  point  to  the 
curve-fitting  analysis.  The  results  are  shown  in  Figures  26,  27,  and  28  for  35, 

98,  and  140  GHz,  respectively. 

At  35  GHz  (Figure  26)  two  things  are  apparent,  h  irst,  the  attenuation  through 
melt-freeze  snow  exceeds  by  a  factor  of  two  or  three  that  through  equitemperature 
snow.  Second,  there  is  larger  scatter  in  the  MF  data.  In  the  case  of  FT  snow,  the 
attenuation  measured  for  vertical  and  horizontal  polarization  have  been  plotted  to¬ 
gether.  This  assumes  polarization  independence  of  attenuation  at  perpendicular 
incidence.  Ai.  average  of  all  the  data  at  each  frequency  bears  this  out. 

For  98  GHz  (Figure  27),  one  sees  an  increase  in  attenuation  over  that  for 
35  GHz.  The  ordinate  scale  has  been  compressed  to  one-half  that  used  for  the 
35-GHz  data.  Again,  the  attenuation  through  MF  snow  is  greater  than  that  through 
FT  snow.  The  parabolic  curve  fit  to  the  melt-freeze  data  is  quite  good  and  shows 
a  leveling-off  as  the  slab  thickness  increases.  On  the  other  hand,  there  is  a  large 
spread  between  the  vertically  and  horizontally  polarized  measurements  through 
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Figure  26.  Attenuation  vs  Slab  Thickness  at  35  GHz, 
(polarizations  vertical  o,  horizontal  x) 
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Figure  27.  Attenuation  vs  Slab  Thickness  at  08  GHz, 
(polarizations  vertical  o,  horizontal  x) 


FT  slabs.  The  curve  fit  to  these  data  is  correspondingly  poorer.  As  the  horizon¬ 
tally  polarized  data  look  better,  a  parabola  has  been  fitted  to  these  points  alone 
(the  dashed  curve  in  Figure  27).  At  140  GHz  (F  igure  28),  the  attenuation  coefficient 
through  MF'  snow  is  again  seen  to  be  greater  than  through  FT  snow. 

Attenuation  coefficients  taken  from  the  slope  of  the  straight  lines  at  35  GHz  and 
from  the  initial  slopes  of  the  parabolas  at  98  and  140  GHz  are  presented  in  Table  8. 


140 GHi  ATTENTION 


MfLT-fREF?t  Eooi’t  MPfPATURE 


THICKNCSSlcml  TmCKWSSIcm) 


I  igure  2o.  Attenuation  vs  Stab  Thickness  .it  140  till/  , 
(polarizations  vertical  o,  horizontal  \l 


At  ;J5  (111/,  the  attenuation  coefficient  measured  through  Ml'  snmi  is  n.i-r.  »ti  .n 
double  that  through  ITT  snow,  there  is  a  sharp  increase  in  I.  as  the  I' requeues  in¬ 
creases  to  98  (ill/.  Again,  the  value  through  Ml  snow  is  greater  than  that  through 
(T1  snow,  the  value  pertaining  to  horizontal  polariz  ition  alone  is  shown  in  pircii- 
tlieses.  At  MO  (ill/,  the  attenuation  coefficient  has  become  even  greater.  At  this 
frequency  it  is  essentially  equal  for  both  types  of  snow. 


Table  li.  Attenuation  Coefficient  vs  Frequence, 
(initial  slope  ol  curves  fitted  to  data) 


1  requeue v 

Attenuation  C 

oefficient  tdli/  1  ■) 

((ill/) 

Alelt-I-  reeze  Snow 

Kquitcmper  dure  Snow 

:ir> 

7  1 

98 

•1  Hi 

l  n:i(:i  l  o' 

MO 

500 

_ 

*95 

if  multiple  scattering  is  much  weaker  at  1  a  (.ill/  (judging  from  the  absent  c  ol 
modulation  on  the  H5-(iHz  data  compared  to  the  higher  frequencies)  tbtm  thicker 
samples  would  be  needed  to  obtain  results  at  this  frequency  comparable  to  the  one 
at  the  higher  frequencies.  To  prove  this,  additional  measurements  are  required. 
Measurements  of  attenuation  should  preferably  be  made  on  a  series  of  snow  slabs 
each  of  increasing  thickness  but  otherwise  unchanging  characteristics. 


4.3  ( loin  par  jgon  Kill)  Ollier  Results 
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In  the  report  discussing  the  in-situ  data,  backscatter  coefficients  in  the 
sample  and  m-situ  modes  are  compared  with  the  results  obtained  by  other  re¬ 
searchers.  In  this  report,  discussion  will  be  limited  to  a  comparison  of  the  slab 
backscatter  data  with  the  in-situ  data.  When  operating  in  the  sample  mode  it  is 
implied  that  cutting  and  removing  the  slab  from  the  snow  cover  does  not  perturb  its 
scattering  and  attenuation  properties.  The  3-dB  beamwidths  of  the  three  trans¬ 
mitting  antennas  are  narrow  enough  so  that  only  a  portion  of  the  slab  is  illuminated. 
Under  this  condition,  the  finite  slab  appears  equal  to  one  infinite  in  size.  Further, 
the  slab  must  be  thick  enough  so  that  the  snow/air  interface  on  its  underside  plays 
no  role  in  the  strength  of  the  scattered  signal.  The  validity  of  these  assumptions 
is  tested  by  comparing  backscatter  coefficients  from  both  types  of  measurements. 
Ideally,  comparisons  should  be  made  on  the  same  day  and  on  the  same  snow.  This 
was  impossible  due  to  the  complexity  of  equipment  reconfiguration  from  one  mode 
to  the  other.  Sample  measurements  were  therefore  carried  out  during  the  earlier 
part  of  the  season,  in-situ  measurements  during  the  later  part.  Sample  data  gathered 

in  February  and  early  March  are  compared  here  with  in-situ  data  gathered  in 
14 

March.  Sample  data  represent  ET  snow  conditions.  These  conditions  persisted 
some  time  into  the  in-situ  measurement  period,  making  the  two  types  more 
comparable.  The  in-situ  program  lasted  through  the  remainder  of  the  winter, 
including  the  late  diurnal  melt-freeze  cycles.  Hence,  some  differences  in  the  re¬ 
sults  may  be  attributable  to  differences  in  the  snow.  In  general,  the  two  modes  of 
operation  produced  similar  output.  Averages  of  <7°  obtained  in  the  two  modes  are 
shown  in  Figures  29,  30,  and  31,  one  figure  for  each  frequency.  The  agreement 
between  the  sample  and  in-situ  measurements  is  best  at  140  GHz,  with  differences 
ranging  to  no  more  than  3  dll.  The  sample  data  show  somewhat  higher  levels  of 
a°VV  than  a°HH'  T  his  is  not  repeated  in  the  in-situ  data  and  we  conclude  that  a 
polarization  dependence  of  the  like-polarized  backscatter  coefficient  cannot  be 
claimed  from  the  available  results.  This  fact  is  corroborated  by  what  is  seen  at 
the  two  lower  frequencies  in  Figures  29  and  30.  Averaging  over  all  like-polarized 
<J°  at  140  GHz  (Figure  31)  leads  to  a  2.  1-dll  difference  between  90  and  45-deg  levels, 
and  to  a  9.  6-dll  difference  between  90-  and  15-deg  levels.  This  matches  with  the 
expected  drop-off  for  incoherent  scatter.  Gross-polarized  backscatter  is  significant 
and  within  3  dB  of  the  like-polarized  values  at  all  angles  for  both  typos  of  measure¬ 
ment. 


Figure  29.  Comparison  of  Sample  and  In-situ  Data  at  35  GHz 


Almost  the  same  good  agreement  is  found  for  98-CHz  results.  There  is  a 
slight  indication  of  a  level  drop  between  90  and  45  deg,  corresponding  to  a  stronger 
specular  component.  The  difference  between  like-  and  cross -polarized  a°  is  also 
larger,  of  the  order  of  4  to  5  dB.  The  tendency  of  the  in-situ  data  to  exceed  the 
sample  data  is  almost  consistent  for  all  types  of  <7°  and  for  all  angles.  Based  on 
previous  reasoning,  and  in  consideration  of  attenuation  coefficients  measured  at 
this  frequency,  we  still  feel  that  differences  must  be  of  a  statistical  nature  as  im¬ 
plied  at  140  GHz  and  not  due  to  the  limited  depth  of  the  slabs. 

Conditions  are  distinctly  different  at  35  GHz.  Here  we  note  a  characteristic 
and  systematic  difference  of  8  to  10  dB  between  the  cross -polarized  backscatter 
coefficient  measured  on  samples  and  on  the  ground.  We  interpret  the  difference 
in  two  ways.  Path  attenuation  through  snow  is  low  enough  so  that  contributions 
from  the  back  side  of  the  slab  definitely  contribute  to  the  cross -polarized  signal. 
Due  to  generally  greater  snow  depth  on  the  ground,  more  scattering  occurs  in  the 
in-situ  case.  There  is  further  the  fact  that  the  ground  below  the  snow  reflects  the 
35 -GHz  signal  which  upon  reflection  illuminates  more  snow  particles  for  more 
cross -polarized  backscatter.  The  like-polarized  data  can  be  used  to  support 
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Figure  30.  Comparison  of  Sample  and  In-situ  Data  at  98  GHz 


this  argument.  At  perpendicular  incidence,  the  snow  surface  is  generally  smooth 
enough  at  35  GHz  to  cause  a  strong  specular  response,  equal  to  or  higher  than  the 
noncoherent  response.  No  substantial  snow  depth  is  required  to  generate  this  sig¬ 
nal  and  correspondingly  a°yy  and  0°^^  arc  comparable  for  the  sample  and  in-situ 
data.  It  is  at  the  lower  grazing  angles  that  the  depth  of  the  snow  comes  to  bear  on 
the  strength  of  the  like-polarized  cr°.  The  sample  data  show  the  characteristic 
specular  response  seen  previously  in  this  report.  That  is.  from  the  very  high 
perpendicular  return  at  90  deg  there  is  an  abrupt  lowering  at  45  deg  to  incoherent 
signal  levels  from  the  slab  which  is  not  thick  enough  to  represent  infinite  snow 
depth  at  this  frequency.  Whether  or  not  an  average  snow  depth  of  30  cm  found 
during  in-situ  measurements  does  this  is  debatable.  There  is  significantly  higher 
like -polari zed  backscatter  than  observed  on  samples.  I  bis  explanation  is  sup¬ 
ported,  if  one  assumes  that  at  lower  grazing  angles  reflection  off  the  ground  becomes 
fairly  effective,  with  Ihe  ground-reflected  signal  again  producing  like -polarized 
backscatter  off  the  snow  cover.  In  summary,  it  appears  to  be  safe  to  conclude  that 
at  98  and  140  GHz  the  sample  measurement  technique  represents  in-situ  conditii  ns 
quite  well.  At  35  GHz,  the  scatter  levels  are  lower  than  what  can  be  expected  lor 
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m-situ  snow  ii  great  depth.  In  nianv  practical  situations  the  snow  depth  is  limited. 
I  inlet-  these  >  onditions  the  slab  data  mav  be  quite  representative  at  this  Irequencv. 

A  review  of  the  subiect  of  nm  rmvave  attenuation  through  snow  cover  is  con¬ 
tained  in  the  report  bv  Stiles  and  Llabv.  lt'  In  addition  to  attenuation  data  thev 
collected  at  d5.  6  GHz.  thev  discuss  attenuation  measurements  made  bv  Currie  et  al 

23 

at  35  and  : 1  5  GHz  and  nieasur-Miients  made  bv  Hatties  and  Crane”  at  do  GHz  using 
snow  produced  in  the  laboratory.  Figure  32,  taken  directly  from  the  I'niversitv- 
ol-Kansas  report  shows  the  procedure  of  snow  attenuation  measurements.  A  com¬ 
posite  layer  is  formed  bv  stalking  slabs  from  the  snow  cover  to  steowise  increase 
the  thickness  of  the  layer.  1  ho  results  of  these  measurements  are  shown  in  f  ig¬ 
ure  33.  1'his  figure  has  also  been  taken  Irom  ihc  Stiles  and  Llabv  report.  Karlier 
it  was  suggested  that  at  do  GHz  our  slabs  were  not  thick  enough  to  show  a  nonlinear 
dependence  of  attenuation  upon  slab  thickness.  This  is  borne  out  bv  the  Stiles  and 
IJlabv  data  for  dry  snow.  A  flattening  in  the  attenuation  curve  in  Figure  dd  appears 


2d.  Battles,  J.W.,  and  Crane,  1).  K.  (11)661  Anemi.it ion  of  Ka -Hand  Knergv  bv 
Snow  and  Ice,  U.  S.  Naval  Ordnance  I  ,ab.  .  NOI.C  Kepi  G7 0. 

Corona,  California,  All  638303. 
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onlv  when  the  slab  thickness  approaches  40  cm.  During  these  measurements  snow 
was  simply  characterized  as  dry,  wet,  or  very  wet.  It  should  be  noted  that  the 
attenuation  coefficient  measured  by  us  through  ET  snow  agrees  fairly  well  with  the 
average  slope  before  the  flattening  of  the  curve  shown  for  dry  snow  in  Figure  33. 
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Figure  32.  Diagram 
Illustrating  the  Stiles 
and  (Jlabyl8  Attenuation 
Measurement  Through 
Lavers  of  Snow  at 
35.  6  GHz 


DATE  J/23-3/?5/TT 


THICKNESS  Of  LAYER  (cm) 


Figure  33.  Path  Loss 
Measured  by  Stiles  and 
L’labv1®  as  a  Function 
of  Snow  Thickness 


Currie's  et  al  experimental  arrangement  consisted  of  two  truck-mounted  pulse 
radars  reflecting  signals  off  a  10-in.  corner  reflector.  The  two-way  attenuation 
through  snow  was  measured  by  placing  slabs  in  front  of  the  corner  reflector  and 
measuring  the  radar  return.  A  reference  signal  was  obtained  bv  raising  the  corner 
reflector  above  the  snow  slab  and  then  measuring  the  return.  Their  attenuation 
data  are  given  in  Table  9.  Snow  properties  in  lheir  report  are  not  described  in  any 
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greater  detail  than  that  listed  in  Table  9.  It  is  therefore  impossible  lo  make  a 
detailed  comparison  of  their  data  and  ours.  On  the  one  occasion  that  thev  made 
simultaneous  measurements  at  35  and  95  GHz,  the  values  of  the  attenuation  coeffi¬ 
cient  thev  obtain  are  in  rough  agreement  with  those  we  measured  through  MF  snow. 
The  attenuation  coefficient  on  15  April  at  35  GHz  seems  anomalous lv  high.  As  ue 
have  seen  in  our  data,  rotation  of  the  snow  slab  results  in  a  random  fluctuation  of 
the  transmitted  signal.  It  is  conceivable  that  this  high  value  is  due  to  critical 
positioning  of  the  slab.  Currie  et  al  measured  an  increase  in  the  attenuation  coeffi¬ 
cient  at  35  GHz  through  wet  snow  but  essentially  no  change  at  95  GHz.  We  did  not 
measure  attenuation  through  wet  snow,  but  as  part  of  our  in-situ  measurements  we 
did  determine  the  effect  of  snow  wetness  upon  the  backscatter  coefficient.  We 
found  a  decrease  in  the  effect  of  free  water  upon  the  backscatter  signal  as  the  fre¬ 
quency  increases. 


Table  9.  Loss  Measurements  at  35  GHz  and  95  GHz  by  Currie  et  al* 


Date 

Time 

Frequency 

(GHz) 

Polariz. 

Snow 

Thickness 
(cm ) 

Snow 

Condition 

Loss 

(dB) 

Atten. 
Coe  I  f. 
(dB/m ) 

14  April  1976 

95 

HH 

6.  0 

Wet 

30 

250 

15  April  1976 

35 

HH 

4.  6 

Frozen 

Snow 

Crust 

21 

228 

16  April  1976 

0635 

35 

HH 

5.  0 

Frozen 

Snow 

Crust 

5 

50 

95 

HH 

5.  0 

Frozen 

Snow 

Crust 

30 

300 

0639 

35 

VV 

5.  0 

10 

100 

07  20 

35 

HH 

3.  5 

3 

43 

07  23 

35 

VV 

3.  5 

7 

100 

0746 

35 

VV 

3.  5 

Melting 

Snow 

Crust 

14 

200 

0748 

35 

HH 

3.  5 

1  1 

157 

Battles  and  Crane  used  an  interferometer  to  measure  attenuation  at  35  GHz 
through  snow  and  ice  manufactured  in  an  environmental  chamber.  The  objective 
of  their  study  was  to  determine  the  effect  on  attenuation  of  the  variation  of  snow 
density  at  a  fixed  temperature.  All  tests  were  made  at  -25°C.  Their  findings  are 
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presented  in  Table  10.  They  claim  their  loss  measurements  are  accurate  to  within 
0.  5  dB.  Their  "snow"  was  produced  by  spraying  a  line  mist  of  cold  water  or  hot 
water  into  the  cooled  environmental  chamber.  By  varying  the  water  temperature, 
pressure,  and  chamber  temperature  they  were  able  to  produce  snow  particles  of 
varying  sizes.  Snow  particles  produced  in  this  way  most  closely  resemble  small 
hail  or  sleet  particles.  The  size  of  the  particles  ranged  from  0.  3  mm  for  the  fine 
snow  up  to  almost  2  mm  for  the  large  ice  crystals.  Samples  composed  of  large  ice 
crystals  are  similar  to  the  Ml-'  slabs  used  in  our  tests.  In  this  case,  the  attenua¬ 
tion  coefficient  (65.7  dB/ml  reported  by  Battles  and  Crane  is  in  rough  agreement 
with  ours  (average  74  dB/ml.  They  also  show  the  effect  of  increasing  the  density 
of  the  snow  by  packing  it.  Scattering  should  decrease  and  consequently  attenuation 
should  decrease.  Also,  if  the  size  of  the  particles  is  small  compared  to  the  wave¬ 
length,  scattering  and  thus  attenuation  should  decrease  likewise.  This  is  true  for 
the  fine,  loose  snow  where  the  ratio  of  particle  diameter  to  wavelength  is  approxi¬ 
mately  0.03.  The  attenuation  coefficient  they  determined  through  ice  is  four  times 
higher  than  the  calculated  value  in  Table  7.  it  is  not  known  if  their  ice  slab  was 
free  from  air  bubbles  which  would  tend  to  increase  attenuation  due  to  scattering 
from  the  bubbles. 


Table  10.  Loss  Caused  by  Snow  and  Ice  at  35  GHz, 
(Battles  and  Crane^) 


Type 

Thickness 

(cm) 

Density 
(g/ cm3) 

— 

Loss 

(dB) 

Atten. 

Coeff. 

(dB/m) 

Fine  loose  snow 

14.  0 

0.  20 

2.  5 

17. 9 

Locally  packed  loose  snow 

14.0 

0.  33 

7.  2 

51.4 

Large  ice  crystals 

14.  0 

0.  39 

9.  2 

65.  7 

Packed  snow 

14.  0 

0.  47 

1.  4 

10.  0 

Ice 

5.  1 

0.  92 

1.  2 

23.  4 

5.  CONCLUSIONS 


Metamorphic  snow  in  the  frozen  state  represents  a  medium  of  low-loss,  dense¬ 
ly  packed  particles  of  resonant  size  in  the  millimeter  wave  region.  Sometimes, 
depending  on  the  process  of  development,  internal  bonds  are  strong  enough  to 
approximate  an  ice  layer  of  relative  smoothness  in  terms  of  wavelength.  It  was  the 
purpose  of  this  report  to  determine  scattering  and  attenuation  characteristics  of  the 
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medium  at  the  atmospheric  -  window  frequencies  ol' 35,  a»,  anti  140  GHz.  Measure¬ 
ments  were  ubtained  on  slabs  of  snow  which  could  be  conveniently  removed  from 
natural  snow  cover.  This  technique  is  almost  mandatory  in  order  to  obtain  attenua¬ 
tion  data,  but  proved  to  be  advantageous  also  in  studying  backscatter  and  bistatic 
scatter  phenomena.  One  major  advantage  of  the  experimental  system  developed 
for  the  program  is  its  ability  to  provide  contiguous  coverage  in  grazing  angle.  Care 
was  taken  to  contain  the  area  of  wave/ snow  interaction  within  the  lateral  extent  of 
the  slab  and  to  keep  the  scatterers  within  the  far  field  of  the  antennas  used.  Snow 
depth  proved  to  be  one  critical  parameter  that  could  not  be  fully  satisfied  at 
35  GHz,  the  lowest  frequency  investigated.  Significant  contributions  from  snow 
depths  greater  than  the  slab  thickness  could  be  expected  as  a  comparison  with  simi¬ 
lar  measurements  on  in-situ  snow  reveals.  Slab  thickness  was  no  problem  at  the 
higher  frequencies,  98  and  140  GHz. 

Two  parameters  were  investigated  regarding  their  dependence  on  grazing  angle 
to  aid  in  the  interpretation  of  the  received  signal  and  to  reveal  the  spatial  arrange¬ 
ment  of  scatterers.  These  were  the  number  of  fluctuations  due  to  phasing  observed 
per  deg  slab  rotation,  and  the  drop-off  in  signal  strength  as  the  grazing  angle  re¬ 
duces.  .Neither  provided  unambiguous  evidence  as  to  the  scatter  mechanisms  in¬ 
volved.  It  was  found  t fiat  reflection  from  a  smooth  dielectric  interlace,  scattering 
from  a  rough  surface,  and  scattering  from  a  volume  of  particles  were  generally 
not  observable  individually.  The  one  exception,  35 -GHz  scatter  from  an  ice -like 
surface  where  specular  reflection  predominates  at  perpendicular  incidence,  is 
somewhat  in  doubt  because  deeper  snow  might  have  resulted  in  a  more  substantial 
volume  scatter  contribution.  Apart  from  diffraction  t  t feels,  which  must  be  con¬ 
sidered  edge  effects  on  the  fluctuation -rate  measurement,  neither  of  the  three 
scatter  mechanisms  could  he  singled  out  clearly  on  the  basis  of  angular  changes  in 
the  number  of  fluctuations  observed.  This  is  not  neeessarilv  a  detriment  but  simply 
an  indication  of  the  complexity  of  the  medium  which  interacts  with  millimeter  waves 
in  a  mixed  mode.  Thus,  it  is  partially  reflecting  in  a  specular  manner  if  snow 
composition,  grazing  angle,  and  frequency  are  appropriate,  and  otherwise,  the 
surface  layer  or  deeper  lavs  are  scattering  incoherently,  again  depending  on  the 
three  parameters  above.  It  must  be  noted  that  incoherent  signal  levels  can  be  quite 
comparable  to  coherent  ones. 

In  order  to  assess  the  influence  of  snow  structure  on  scattering  and  attenuation, 
three  types  of  metamorphic  snow  were  identified,  related  to  time  periods  of  the 
progressing  winter  season.  Distinctions  were  made  primarily  in  line  with  meteor¬ 
ological  classifications.  The  measured  data  suggest  some  systematic  differences 
in  millimeter  wave  interaction  with  the  three  snow  types  specified.  Tor  a  svsfem 
designer  they  might  be  largely  unimportant.  A  convincing  difference  between 
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vertical  and  horizontal -polarization  results  could  not  be  established  at  any  fre¬ 
quency  or  any  grazing  angle  investigated.  This  applies  to  backscatter  and  attenua¬ 
tion.  Cross -polarized  backscatter  was  surprisingly  high,  reaching  levels  within 
several  dB  of  like -polarized  backscatter,  when  the  incoherent  scatter  mechanism 
prevailed.  When  specular  scatter  dominated,  the  cross -polarized  component  was 
low. 

At  90  deg  grazing  angle,  average  backscatter  coefficients  were  measured  at 
35,  98,  and  140GHz:  off  sleet  11,  9.  11.5  dB,  off  MF  snow  3,  8,  9.  5  dB,  off 
ET  snow  6,  4.5,  9  dB,  off  ET  snow  (cross-polarized)  -13,  -1,  7.5  dB.  At  15  deg 
grazing  angle,  average  backscatter  coefficients  were  measured  at  35,  98,  and  140  GHz 
off  sleet  -14,  -6,  -3  dB,  offMF  snow  -11,  -1.5,  -0.  5  dB,  offET  snow  - 13,  -4.5, 

0.5  dB,  off  ET  snow  (cross-polarized)  -24,  -9.5,  -3  dB.  The  attenuation  coefficient 
at  90  deg  penetration  into  the  slab  measured  at  35,  98,  and  140  GHz:  through 
MF  snow  74,  413,  and  506  dB/m  and  through  ET  snow  29,  183  (318)  and  495  dB/m. 
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Appendix  A 

Snow  Data 


Physical  parameters  of  the  snow  cover,  from  which  snow  slabs  were  cut,  are 
presented  in  this  appendix.  One  combination  table/diagram  is  given  for  every  day 
on  which  millimeter-wave  measurements  were  conducted  (Tables  A1  through  A20). 
The  data  characterize  the  snow  at  the  time  the  slabs  were  cut.  A  standard  form  of 
presentation  has  been  used.  The  vertical  scale  at  the  center  of  the  form  gives  the 
height  above  ground.  Any  snow  parameter  indicated  on  the  form  at  this  height  is 
pt.taining  to  this  particular  height.  The  stratigraphy  of  the  snow  cover  is  shown  by 
hor.zontal  lines  drawn  between  the  left  edge  of  the  form  and  the  height  scale.  Most 
of  the  column  headings  are  self-explanatory.  The  "I.aver"  column  normally  con¬ 
tains  a  one-word  description  of  the  layer  such  as  hard,  icv,  brittle,  and  so  on.  In 
the  "Density"  column,  a  tick  mark  is  placed  to  show’  the  height  at  which  the  measure¬ 
ment  was  made.  The  air  temperature,  indicated  by  the  *  symbol  on  the  tempera¬ 
ture  plot,  was  obtained  at  a  position  just  above  the  surface  of  the  snow.  Under 
"Comments"  is  placed  such  information  as  surface  hardness,  weather  data,  data  on 
the  history  of  the  snow,  and  in  some  cases  additional  data  on  the  stratigraphy  of 
the  snow.  The  snow  was  completely  frozen  at  the  time  the  test  slabs  were  cut.  No 
melting  occurred  during  the  tests.  Free-water  content  is  therefore  not  listed  on 
the  form. 
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Table  A2.  Physical  Parameters  of  Snow  Cover  17  January  1978 


Table  A3.  Physical  Parameters  of  Snow  Cover  19  January  1978 


Table  A4.  Physical  Parameters  of  Snow  Cover  23  January  1978 


Hounded 


Table  A7.  Physical  Parameters  of  Snow  Cover  28  January  197 


able  A8.  Physical  Parameters  of  Snow  Cover  30  January  1978 
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Table  A9.  Physical  Parameters  of  Snow  Cover  31  January  1978 


Table  A10.  Physical  Parameters  of  Snow  Cover  1  February  1978 
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Sky  was  overcast  during  most 
of  day.  Air  temperature  did 
not  rise  above  -6°C. 


Table  A12.  Physical  Parameters  of  Snow  Cover  4  February  1978 


Table  A13.  Physical  Parameters  of  Snow  Cover  21  February  1978 
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r  a  meters  of  bnow  Cove 


ea  rlier 


Physical  Parameters  of  Snow  Cover  24  February  197  8 


ai.d 

earlie  r 
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Table  A20.  Physical  Parameters  of  Snow  Cover  2  March  1978 


earlier 
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